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SECTION  I 


INTRODUCTION 

The  Ground  Electro-Optical  Oeep  Space  Surveillance  System 
(GFODSS)  is  tasked  to  detect  faint  space  objects  in  orbits  exceeding 
3000  miles  altitude.  The  surveillance  function  is  performed  by 
optical  telescopes  equipped  with  electro-on tical  sensors.  In  this 
application  the  sensor  has  to  combine  high  sensitivity  with  high 
resolution  and  a large  number  of  independent  resolution  cells. 

All  these  requirements  call  for  a TV  type  sensor;  the  specific 
device  chosen  for  the  GEODSS  system  was  the  FBSICON  image  tube. 

The  FBSICON  Sensor  tube  (Electron  Bombarded  Silicon  Target 
Sensor  Tube)  combines  high  sensitivity  with  high  resolution^  while 
maintaining  mechanical  structural  compactness  and  simplicity.  In 
the  special  case  of  space  surveillance  systems,  the  high  resolution 
allows  for  the  observation  of  faint  ob-jects  in  large  fields,  while 
the  high  sensitivity  permits  operation  in  the  photon-no-*  se  limited 
domain. 

Television  systems  in  general  are  designed  to  observe  scenes 
with  limited  ranges  in  illumination  and  contrast.  Many  of  the 
operational  parameters  of  a given  TV  system  are  specifically  chosen 
to  maximize  the  output  for  human  observers  watching  the  scenes  on 
phosphor  screens.  The  performance  specifications  furnished  by  the 
manufacturers  consist  mainly  of  signal’-trans fer  and  modulation  transfer 
curves.  This  information  is  usually  quite  sufficient  to  design  com- 
mercial TV  systems  and,  therefore,  has  become  standardized. 
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In  astronomy  and  space  surveillance,  however,  the  objects  of 
interest  are  point  images  of  very  different  intensities,  riding  on  a 
relatively  uniform  background.  In  many  cases  the  video  signal  is  pro- 
cessed directly  and,  therefore,  scan  patterns,  integration  times  and 
readout  times  can  be  optimized  for  any  given  operation.  The 

specific  application  of  EBSICON  sensor  tubes  in  the  GEODSS  system 
requires  the  knowledge  of  tube  performance  data  for  a wide  variety 
of  operational  parameters.  Furthermore,  the  performance  of  sensor 
tubes  not  yet  in  existence  must  be  projected  and  their  usefulness 
in  various  system  configurations  must  be  analyzed. 

The  signal-transfer  and  sensor  modulation-transfer  functions 
do  not  adequately  measure  the  sensor  tube  performance  in  this  appli- 
cation or  the  response  characteristics  for  point  images.  Furthermore, 
the  literature  does  not  provide  guidelines  for  extrapolating  per- 
formance characteristics  from  one  set  of  given  operational  parameters 
to  another.  To  accomplish  the  tasks  specifically  required  in  the 
GEODSS  program,  it  was  necessary  to  develop  a mathematical  model  as 
a tool,  based  on  the  physical  effects  taking  place  inside  the  sensor; 
and  to  trace  the  signal  and  the  noise  through  the  different  stages 
of  the  sensor,  up  to  and  including  the  video  signal  output.  This 
mathematical  model  can  be  used  to  aid  the  laboratory  test  programs  and 
tube  development  programs  by  providing  a proper  theoretical  basis  for 
explanation  of  observed  effects  and  by  predicting  the  results  of 
the  experiments. 
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The  tube  model,  as  it  is  presented  in  this  report,  provides  signal 
transfer  characteristics  both  for  point  sources  and  flat  fields.  It 
also  provides  the  noise  transfer  characteristics  of  EBSICON  tubes  for 
photon  noise  limited  inputs  for  a large  variety  of  operational 
parameters.  The  noise  characteristics  of  the  output  are  described 
in  terms  of  the  actual  density  distribution  functions  both  for  the 
signal  and  noise.  These  are  of  special  importance  for  calculation 
of  detection  probability  and  false  alarms.  The  model  gives  the 
proper  lag  data  and  such  special  effects  as  signal  splitting  in 
point  images.  The  model  is  quantitative,  i.e.,  all  the  physical 
parameters,  such  as  material  constants  and  physical  dimensions,  are 
chosen  to  very  closely  represent  actual  sensor  tubes. 

This  report  consists  of  four  major  parts.  The  mechanical 
layout  and  electrical  operation  is  described  to  put  the  mathematical 
model  in  perspective.  The  physical  effects  taking  place  in  the 
sensor  are  analyzed  to  a level  of  detail  which  is  required  to 
delineate  the  mathematical  model.  The  physical  effects  are  trans- 
lated into  the  mathematical  model;  the  equations  derived  and  solved. 

And  finally,  the  results,  as  they  are  obtained  by  the  mathematical 
model,  are  presented  and  compared  with  laboratory  results. 
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SECTION  II 


DESCRIPTION  OF  SENSOR  TUBE 

2.1  MECHANICAL  LAYOUT 

2.1.1  General 

An  EBSICON  tube  consists  basically  of  three  parts,  all  enclosed 
in  one  single  vacuum  envelope:  The  image  section,  the  storage  target 

and  the  read  section.  The  overall  dimensions  of  the  envelope  vary, 

( ependino-  on  the  specific  tube  model  and  manufacturer.  A separate 
cylindrical  magnetic  focusing  and  deflection  coil  assembly  surrounds 
the  read  section.  A detailed  outline  and  a photograph  of  a typical 
EBSICON  are  shown  in  Figures  II-l  and  II-2  respectively.  Except  for 
the  field  mesh  and  the  heater  of  the  read  gun,  the  tubes  have  complete 
rotational  symmetry.  All  electrodes  of  the  image  section  and  the 
storage  target  assembly  are  mounted  on  circular  disks  which  are  held 
in  place  either  by  glass  or  ceramic  spacer  rings.  These  spacer  rings 
are  braced  to  the  circular  mounting  disks  of  the  electrodes,  and  form 
the  actual  vacuum  envelope.  The  mounting  disks  protrude  beyond 
the  spacers  in  order  to  apply  the  necessary  voltages  directly.  The 
electron  gun  is  a compact  assembly  of  tightly  spaced  electrodes 
connected  to  a standard  pin-type  tube  socket,  which  is  inserted  as 
a subunit  into  the  end  of  the  read  section.  This  type  of  design 
provides  for  a very  rugged  and  clean  device,  with  short  electrical 
connections  and  no  mechanical  adjustments,  except  for  the  placement  of 
the  magnetic  coil  assembly. 
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Figure  U- 1 OUTLINE  OF  TYPICAL  EBSICON 
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Figure  II-2  PHOTOGRAPH  OF  TYPICAL  EBSICON 


2.1.2  The  Image  Section 


The  image  section  has  the  function  of  converting  the  spatially- 
distributed  optical  image  at  the  photocathode  into  a charge  image 
which  is  stored  in  the  target  and  subsequently  read  out  electronically. 

The  optical  image  on  the  photocathode  causes  photo  electrons  to  be 
released  from  the  photocathode.  The  number  of  photoelectrons  released 
depend  on  the  intensity  of  the  optical  radiation.  The  photo- sens i tive 
material  of  the  photocathode  - typically  an  S20  or  S25  - is  deposited 
on  the  inside  of  the  window  of  the  vacuum  envelope  in  order  to  allow 
the  released  photoelectrons  to  enter  the  vacuum  section  directly  and 
to  be  imaged  on  the  storage  target.  The  imaging  function  is  performed 
by  the  focusing  electrode  placed  opposite  the  cathode,  a cone  or 
funnel-shaped  piece  of  metal,  which  defines  the  electrostatic  focusing 
field.  Some  image  sections  are  equipped  with  additional  electrodes 
to  permit  zooming  or  gating  of  the  image.  Behind  the  focusing  electrodes 
is  the  storage  target  assembly,  which  holds  the  storage  target  and 
some  additional  components  such  as  the  target  mesh. 

A fiber-optic  faceplate  is  attached  to  the  vacuum  envelope  and 
serves  as  the  entrance  window,  providing  a flat  surface  for  the  optical 
image-plane  interface  and  a curved  surface  on  the  other  side  to  hold 
the  photocathode.  This  curvature  is  necessary  to  correct  for 
aberration  in  the  image  section,  which  in  the  case  of  a flat  cathode  would 
reduce  the  resolution  significantly. 
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2.1.3  Hie  Storage  Target 


The  storage  target  of  an  EBSICON  tube  is  a thin  silicon  wafer  which 
carries  an  array  of  small  diodes.  The  detailed  target  structure  and 
the  diode  spacing  differs  depending  on  the  manufacturer;  in  this  report 
the  most  advanced  target  design,  the  Westinghouse  ’’Deep  Etched  Go  Id 
Capped  Diode  Array  Target"  will  be  described.  This  type  of  target 
structure  originated  with  Westinghouse  Corporation  and  was  chosen  to 
become  the  target  for  the  GFODSS  sensor.  Sponsored  by  the  Electronic 
Systems  Division  of  the  Air  Force,  with  the  technical  assistance  of  the 
Army  Night  Vision  Laboratory,  an  EBSICON  sensor  tube  development 
program  was  initiated*  Texas  Instruments  is  performing  the  target 
development  and  Westinghouse  Corporation  is  handling  the  tube  design 
to  incorporate  this  type  of  target. 

A square  raster  of  deep  furrows  is  etched  into  the  silicon  wafer, 
which  leaves  a number  of  square  "mesas".  These  mesas  are  doped  to 
form  n-p  diodes  and  have  a solid  gold  cap  attached,  which  covers  the 
diode  itself  and  most  of  the  space  between  the  diodes.  The  diode 
spacing  is  12.5  pm  center-to-center  (2000  diodes/inch);  the  gap 
between  gold  caps  is  approximately  2-3  ym.  A reflecting  n+  layer 
and  a separate  buffer  layer  are  added  to  the  bulk  of  the  silicon 
wafer  on  the  opposite  side  from  the  diodes.  Schematic  diagrams  of 
the  target  structure  and  of  the  diode  array  are  shown  in  Figure  II-3 
and  Figure  II-4  respectively. 
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Figure  IT- 3 SCHEMATIC  DIAGRAM  OF  TARGET  STRUCTURE 
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Figure  n-4  SCHEMATIC  OF  DIODE  ARRAY 
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2.1.4  The  Read  Section 


The  read  section  consists  of  the  electron  gun,  the  electron 
drift  space  and  the  magnetic  deflection  coils.  The  electron  gun 
resembles  a triode  in  principle;  electrons  emitted  by  a thermal 
cathode  are  focussed  by  a negative  control  electrode  G^  and  acceler- 
ated by  a plate  at  positive  potential,  . A fraction  of  the  electrons 
passes  through  a small  aperture  in  the  positive  plate  and  forms  a well- 
defined,  very  narrow  electron  beam.  Typical  parameters  of  the  beam 
are  approximately  50  pm  diameter  and  1 to  3 y amp  beam  current. 

This  electron  beam  enters  the  drift  space.  In  this  drift  space  the 
magnetic  deflection  of  the  beam  takes  place;  the  beam  can  be  deflected 
both  horizontally  and  vertically  according  to  the  applied  magnetic 
field.  The  field  mesh  at  the  end  of  the  drift  space  defines  an 
equipotential  surface  opposite  the  target. 

In  this  space  between  the  field  mesh  and  the  storage  target,  the 
electrons  are  decelerated  and  land  on  the  target  with  velocities 
close  to  zero.  Since  the  beam  scans  the  whole  target,  the  spatial 
uniformity  of  this  deceleration  field  between  target  and  field  mesh 
has  to  be  maintained  with  high  accuracy.  A schematic  drawing  of  the 
read  section  is  presented  in  Figure  II-5. 

2.2  ELECTRICAL  OPERATION 
2.2.1  General 

In  this  section  the  general  principles  for  operating  an  image 
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Figure  IE- 5 BASIC  STRUCTURE  OF  READ  SECTION 


tube  will  be  discussed.  A detailed  description  of  the  physical 
effects  taking  place  at  the  different  components  of  the  tube  will 
be  presented  in  Section  III  of  this  report.  Here  they  will  be 
sketched  only  briefly,  in  order  to  aid  in  the  understanding  of  the 
general  sensor  operation. 

2.2.2  Electrical  Potentials 

The  various  electrical  potentials  as  they  are  applied  to  the 
electrodes  of  the  sensor  tube  are  listed  below: 


Electro  de: 


Thermal  cathode 


Function: 


Emits  electrons  of 
read  beam 


G 


1 


Focuses  and  controls 
read  beam 


G 


2 


Acelerates  read  beam 
electrons 


Potential: 
Reference  potential 

-50  to  -150^ 

+300  to+350V 


G 


3 


G 


4 


Storage  target 


Photocathode 


Defines  potential  in  +750  to  800V 

drift  space 


Defines  uniformity  of 
field  gradient  between 
drift  space  and  storage 
target 

Stores  photoelectrons 
for  readout 


+800  to  850V 


+10V  bias  on 
illumination  side 
0 to  +10V  on  read 
side 


Emits  photoelectrons  -14  to  -10KV 

depending  on  illumination 
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A potential  profile  along  the  axis  of  the  sensor  tube  is  given  in 
Figure  II-6 . 

The  following  paragraphs  give  a brief  account  of  the  actual 
sensor  tube  functions  in  order  to  present  a basis  for  the  discussion 
of  the  physical  effects  analyzed  in  the  next  section. 

The  incident  photon  flux  releases  photoelectrons  from  the 
photocathode  (which  is  at  4,000  to  10,000  volts  negative).  The 
electrons  are  accelerated  and  hit  the  target  with  high  velocity , 
releasing  additional  charge  carriers  in  the  bulk  of  the  storage 
target.  The  storage  target  itself  is  basically  a capacitor 
whose  surface  facing  the  photocathode  is  connected  to  a low  positive 
DC  bias  (7-10  volt) . The  other  surface  of  the  target  (diode 
side)  is  left  floating  . The  first  scan  of  the  read  beam  puts 
this  surface  at  the  potential  of  the  gun  cathode.  This  results  in 
a voltage  differential  across  the  target  equal  to  the  applied  DC 
bias.  If  photoelectrons  reach  the  target,  the  released  charge 
carriers  discharge  the  target  capacitors,  therefore  driving  the  diode 
surface  which  faces  the  read  beam  towards  the  level  of  the  positive 
DC  bias.  During  the  following  scan,  electrons  from  the  read  beam 
find  a positively  charged  surface  and  land  on  it  until  this  surface 
has  again  reached  gun  cathode  potential.  The  resulting  current  is 
the  observed  signal.  Since  the  number  of  emitted  photoelectrons 
the  number  of  released  charge  carriers,  the  reduction  in  the  voltage 
differential  across  the  storage  target, and  the  number  of  electrons 
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Figure  31  - 6 TYPICAL  POTENTIAL  DISTRIBUTION  ALONG  AXIS  OF  EBSICON  TUBE 


recharging  the  storage  target  are  all  proportional  to  each  other, 
the  resulting  signal  current  is  a linear  function  of  the  incident 
light  intensity.  In  a rough  analogy,  the  storage  target  can  be 
considered  to  be  equivalent  to  a large  number  of  individual  capa- 
citors in  parallel,  each  capacitor  having  a varible,  illumination- 
dependent  shunt  resistor  with  the  read  beam  acting  as  a multiplexer 
switch. 

2.2.3  The  Scan  Process 

The  electron  beam  scans  the  storage  target  in  raster  fashion. 
During  readout  the  read  beam  travels  across  the  storage  target  in  a 
horizontal  line.  After  each  line  scan  the  beam  is  retraced  very 
fast  compared  to  the  horizontal  scan  velocity  and  at  the  same  time 
is  deflected  downward  to  start  a new  horizontal  line.  During  a given 
time,  At,  usually  referred  to  as  the  frame  time,  the  whole  target 
area  is  scanned  and  the  read  beam  is  returned  to  the  starting  point 
of  the  raster.  The  number  of  scan  lines,  the  size  of  the  scanned 
target  area^  and  the  frame  can  be  varied  over  a wide  range  of  values. 

In  commercial  television  systems,  the  sensor  tube  is  operated 
in  the  continuous  mode,  i.e. , the  target  is  exposed  to  photoelectrons 
continuously  during  the  scanning  process.  In  this  case  the  effective 
integration  time  A T is  equal  to  the  frame  time  Ax.  For  astronomical 
and  space  surveillance  work  the  snapshot  or  integrating  mode  is  of 
more  interest.  In  this  case  the  target  is  exposed  to  photoelectrons 
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from  the  photocathode  for  an  exposure  time  AT  (integration  time)  while 
the  scanning  electron  beam  is  turned  off.  After  the  target  is 
exposed,  the  target  is  read  out  during  a read  or  frame  time  At.  This 
mode  of  operation  is  usually  required  where  the  incident  light 
intensities  are  too  low  to  produce  reasonable  signal  levels  in  the 

continuous  mode.  The  ratio  AIT  can  be  set  arbitrarily. 

At 

2.2.4  The  Composite  Video  Signal 

The  scanning  process  converts  the  time- varying,  two-dimensional 
energy  distribution  into  a time-varying  current.  In  order  to  recon- 
struct the  original  image,  the  individual  scan  lines  have  to  be 
indexed  accurately.  After  the  signal  leaves  the  preamplifier 
(which  is  mainly  used  to  match  impedances) , timing  pulses  of  opposite 
polarity  are  added  to  the  signal,  as  well  as  a pulse  to  block  any 
signal  during  line  retrace.  The  pulses  are  adjusted  to  define  a 
black  level,  from  which  the  signal  currents  or  voltages  are  measured. 
An  example  of  a typical  video  output  signal  is  given  in  Figure  II-7. 

For  the  purpose  of  this  report,  it  is  always  assumed  that  the 
electrical  bandwidth  of  the  amplifier  chain  is  sufficient  to  fully 
maintain  the  optical  resolution  of  the  device.  No  degradation  of 
the  optical  resolution  due  to  other  limitations  than  that  of  the 
device  itself  are  considered  here. 
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Figure  IL-7  TYPICAL  STRUCTURE  OF  VIDEO  SIGNAL 


SECTION  III 


PHYSICAL  PHENOMENA  IN  IMAGE  TUBES 

In  order  to  develop  a credible  mathematical  model  for  the  sensor 
tube  performance,  it  is  necessary  to  analyze  the  actual  physical 
effects  taking  place  at  various  sections  of  the  sensor.  This  is 
especially  important  if  a detailed  statistical  treatment  of  the  signal 
transfer  is  required.  The  gain,  storage  and  scan  processes  deter- 
mine the  statistical  distributions  of  electrons  and  have  to  be  incor- 
porated into  a final  noise  analysis.  For  this  reason,  four  main 
areas  are  identified,  where  major  changes  in  electron  distributions 
occur.  These  are  the  photocathode,  the  gain  effects  in  the  target, 
the  electron  storage  process  and  the  read  process.  The  physical 
effects  taking  place  in  these  areas  are  discussed  in  the  following 
paragraphs . 

3.1  THE  PHOTOCATHODE 

The  optical  imaging  system  forms  an  image  on  the  faceplate  of 
the  sensor  tube.  The  optical  image  is  transfered  by  the  fiber- 
optics  directly  to  the  photocathode.  There  the  optical  flux  (photons) 
causes  the  release  of  photoelectrons  in  numbers  proportional  to  the 
intensity  of  the  incoming  radiation.  A certain  number  of  the 
photoelectrons  physically  leave  the  bulk  of  the  photocathode  and 
reach  the  vacuum  of  the  image  section.  The  photon-photoelectron 
transfer  takes  place  extremely  fast;  the  estimated  times  are  less 
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than  1 x 10  sec.  Nonetheless  losses  occur  in  the  process;  they 
can  be  classified  as  photon  and  electron  losses.  A certain  number  of 
the  incident  photons  can  penetrate  the  photocathode  without  any  effect, 
or  are  reflected  by  the  surface  of  the  photocathode.  Another  fraction 
of  the  total  number  of  photons  is  absorbed  in  the  bulk  of  the  photo- 
cathode, but  fail  to  release  photoelectrons.  On  the  other  hand, 
photoelectrons  created  too  far  in  the  bulk  of  the  photocathode  recom- 
bine before  they  can  reach  the  surface  of  the  photocathode. 

Because  of  these  effects,  the  efficiency  of  the  photon-photo- 
electron conversion  process  is  considerably  less  than  unity.  This 
ratio,  the  number  of  released  photoelectrons  to  the  total  number  of 
incoming  photons,  is  called  the  "quantum  efficiency".  This  quantum 
efficiency  is  strongly  dependent  on  the  energy  of  the  incoming 
photons  (spectral  range  of  the  radiation),  but  for  all  practical 
purposes  an  average  quantum  efficiency  can  be  defined  and  used  for 
the  visible  part  of  the  optical  spectrum.  For  photocathodes 
employed  in  modem  sensors  (S20  and  S25)  this  average  quantum 
efficiency  varies  between  5%  and  20%.  This  implies  that  a large 
number  of  photons  is  lost  in  the  conversation  process  and  that  in  an 
analysis  of  the  statistical  distribution  of  the  photoelectron 
counts,  these  loss  effects  can  not  be  neglected. 

The  only  physically  observable  quantities  connected  with  the 
photoelectric  effect  are  the  mean  of  the  photon  flux,  the  mean  of  the 
photoelectron  flux, and  the  distribution  of  the  photoelectron 
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counts  during  any  diven  integration  time.  The  statistical  dis- 
tribution of  the  photons  cannot  be  observed  directly. 

The  statistical  distribution  of  the  photons  is  transferred 
with  modifications  into  the  statistics  of  the  photoelectron  counts. 

The  dominant  factor  in  determining  the  statistics  of  the 
electron  counts  is  the  ratio  of  the  sensor  integration  time  to  the 
coherence  time  of  the  incident  photons.  A detailed  discussion  of 
these  statistics  can  be  found  in  References  1 and  2.  Three  different 
domains  can  be  defined:  1)  If  the  integration  time  (observation 

time)  is  large  compared  to  the  coherence  time,  the  observed 
statistical  distributions  in  the  photoelectron  count  follow  Poisson 
statistics.  2)  If  the  observation  time  is  approximately  equal  to 
the  coherence  time,  the  classical  Bose  - Einstein  distribution  will 
be  observed.  3)  If  the  observation  time  is  less  than  the  coherence 
time  of  the  incoming  photons,  then  a complex  quantum-mechanical 
description  is  required.  Electro-optical  photon-counting  sensors > 
such  as  the  EBSICON,  always  have  large  integration  times  compared 
to  the  coherence  time  and  therefore  the  resultant  observed  photo- 
electron distribution  will  be  Poisson. 

3.2  THE  GAIN  MECHANISM 

In  the  image  section,  the  photoelectrons  are  accelerated  up  to 
10  KeV  and  hit  the  n-type  silicon  of  the  storage  target.  The  high- 
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energy  electrons  penetrate  several  micrometers  and  come  to  rest. 

The  interaction  of  the  high-energy  electrons  with  the  lattice  results 
mainly  in  the  creation  of  electron-hole  pairs  in  the  crystal,  which 
has  a well-defined  threshold  of  3.5  eV.  Beyond  this  threshold  the 
electron  hole-pair  combinations  cannot  carry  any  significant  amount 
of  energy,  therefore  the  energy  dissipation  of  the  high-energy 
photoeTectrons  occurs  in  a very  orderly,  highly  controlled  way: 

In  a step-by-step  procedure  each  creation  of  an  electron-hole  pair 
subtracts  3.5  eV  from  the  energy  of  the  photoelectrons,  until  the 
energy  of  the  photoelectrons  is  exhausted.  The  number  of  electron- 
hole  pairs  per  high-energy  photoelectron  is  called  the  primary  gain 
and  is  the  ratio 


E 

%=t 


(3.2.1) 


where  E is  the  energy  of  the  high-energy  photoelectrons  and  E 
P s 

the  energy  of  the  secondary  electron-hole  pairs.  For  E^  - 10,000  eV 

and  E = 3.5  eV  the  primary  gain  becomes  approximately  3000.  Losses 
s 

occur  due  to  absorption  effects.  A thin  aluminum  layer  is  deposited 
on  the  illumination  side  of  the  target.  This  layer  acts  as  a light 
shield  and  prevents  photons  which  have  passed  through  the  photocathode 
from  entering  the  target  bulk  material.  High-energy  photoelectrons 
penetrate  this  aluminum  layer,  but  lose  some  of  their  energy.  If 
the  layer  is  thin  enough  and  uniform,  the  only  losses  occur  in  the 
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energy  of  the  photoelectrons;  the  number  stays  constant.  When  the 
electrons  reach  the  silicon  bulk  material,  their  effective  energy 
is  that  corresponding  to  the  acceleration  voltage  in  the  image 
section  but  reduced  by  an  energy  loss,  E . This  has  a significant 

Li 

effect  on  the  primary  gain,  G^,  and  equation  (3.2.1)  can  now  be 
rewritten 


G 

P 


E -E 

-P  L 


E 

s 


(3.2.2) 


The  other  important  loss  effects  in  the  target  are  mainly  recombina- 
tions of  charge  carriers  before  they  reach  the  depletion  layer  in 
the  p-n  junction.  The  total  number  of  secondary  charge  carriers 
is  reduced,  and  therefore  also  the  total  gain  of  the  target.  If  one 
defines  the  effective  gain  G^  as  the  primary  gain  G^  reduced  by 
the  loss  factor  L 


G = G (1-L) 
e p 


(3.2.3) 


one  obtains  from  equation  (3.2.2) 

E — E 

G = ^ (1-L)  (3.2.4) 

e E 

s 

The  changes  in  the  statistical  distribution  of  the  charge  carriers 
due  to  the  gain  process  and  the  associated  loss  mechanism  are 
significant.  The  electron  hole-pair  creation  in  the  silicon  bulk 
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constitutes  a highly  deterministic  process  with  insignificant 
statistical  scatter  of  the  number  of  electron-hole  pairs  for  any 
given  photoelectron  energy*  This  in  conjunction  with  uniformity  of 
the  acceleration  field  in  the  image  section,  which  is  required  for  good 
optical  resolution,  guarantees  the  uniformity  of  the  photoelectron 
energies;  therefore  the  statistical  scatter  in  the  total  nuirber  of 
electron-hole  pairs  is  also  insignificant*  The  original  fluctuations 
in  the  number  of  photoelectrons  released  by  the  photocathode  are 
maintained;  their  random  arrival  corresponds  to  packets  of  electron- 
hole  pairs  each.  Superimposed  on  this  distribution  are  the  distributions 
of  the  loss  mechanisms.  The  effects  leading  to  the  losses  E^  (electron 
penetration  of  the  aluminum  light  shield)  are  basically  the  same,  which 
create  electron-hole  pairs  in  the  silicon:  Energy  is  dissipated  in  small 

and  uniform  quantities,  again  a highly  deterministic  process,  which 
does  not  change  the  original  distribution  in  any  significant  way*  Finally, 
the  charge  carrier  recombination  effects  are  a random  process  which,  for 
lossy  targets,  can  result  in  significant  signal  loss*  The  recombination 
losses  for  EBSICON  targets  are  small  and  the  corresponding  statistical 
process  is  essentially  deterministic. 

A detailed  mathematical  treatment  of  the  noise  problems  associated 
with  the  gain  mechanism  is  presented  in  Section  4.1. 

3*3  PHYSICAL  EFFECTS  CONNECTED  WITH  CHARGE  STORAGE 

The  main  function  of  the  storage  target  is  the  integration 
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of  the  charge  created  by  the  photoelectrons  and  to  store  this  charge 
until  it  is  read  out  by  the  read  beam*  The  electron-hole  pairs 
created  by  the  ph'otoelectrons  in  the  n-type  silicon  become  separ- 
ated. The  electrons  move  towards  the  n+  layer  and  play  no  further 
role  in  the  process.  The  positive  holes  (p-charges)  are  extracted 
by  the  depletion  layer  of  the  p-n  junction.  The  depletion  layer 
has  an  intrinsic  capacitance  with  very  low  conductivity,  which 
enables  this  layer  to  maintain  a potential  differential  over  an 
extended  period  of  time.  The  storage  target  represents  basically 
a large  array  of  individual  capacitors  each  originally  charged  with 
a voltage  equivalent  to  the  target  bias  voltage  V_. 

D 

The  positive  charge  carriers  partially  discharge  this 

potential  to  a new  potential  V . The  potential  difference  V - V 

I j B L 

is  a function  of  the  depletion  layer  capacitance  C(V)  and 
the  number  of  charge  carriers.  If  the  total  charge  of  all  charge 
carriers  collected  on  one  diode  is  the  potential  difference 
becomes 


V = v -V  » — 
C B L C 


(3.3.1) 


V 


c 


The  spatial  distribution  of  the  optical  signal  is  now  trans- 
ferred into  voltage  differences  of  a large,  two-dimensional  array 
of  capacitors,  where  the  largest  potential  difference  V^,  corresponds 

to  no  optical  signal  (black)  and  zero  potential  difference  to 
maximum  optical  signal  (white) . 
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The  width  of  the  depletion  layer  is  basically  proportional  to 
the  square  root  of  the  external  voltage  and  since  the  capacitance  C 
of  the  depletion  layer  is  inversely  proportional  to  the  width  of 
the  depletion  layer,  the  change  of  capacitance  C with  voltage  can 
be  expressed  by 


C(V)  = c 

o 


y i + v/v 


(3.3.2) 


where  is  the  capacitance  without  external  voltage,  V,  the  exter- 
nal voltage  and  v,  the  contact  potential.  For  a single  diode  of 

-14 

a deep-etched  target  the  values  for  C are  2 x 10  F and 

o 

0.4  volt.  For  normal  operation  the  changes  in  capacitance  are  not 
critical.  For  voltages  large  compared  with  the  contact  potential 
the’ change  in  capacitance  is  proportional  to  the  square  root  of  the 
voltages,  i.e.,  over  a range  from  10  volt  to  3 volt  the  capa- 
citance increases  only  by  50%.  Furthermore  the  most  rapid  changes 
in  capacitance  take  place  when  the  capacitors  are  practically  dis- 
charged, that  is  at  very  high  light  levels  (Vc  ^v).  In  this  case 
the  nonlinearities  between  the  stored  charge  and  the  resulting 
voltages  (equation  3.3.2)  are  usually  masked  by  other  effects.  For 
these  reasons  it  is  valid  to  assume  an  average  capacitance  for  the 
depletion  layer  and  neglect  its  changes  with  voltage. 
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It  is  obvious  that  the  capacitors  can  only  be  discharged  to 
zero,  where  the  depletion  layer  has  its  smallest  width.  In  case 
of  overload,  i.e.,  the  incoming  optical  signal  had  created  much  more 
positive  charge  carriers  than  are  necessary  for  discharging  the 
capacitor,  the  charge  carriers  are  attracted  by  the  next  diodes 
surrounding  the  discharged  one.  The  charge  carriers  discharge  these 
diodes  in  turn,  and  this  effect  repeats  itself  until  all  charge 
carriers  have  either  reached  a diode  or  have  recombined  in  the  bulk 
of  the  n-layer.  This  effect  is  called  ’’blooming"  which  is  particularly 
noticeable  in  the  case  where  very  bright  point  images  are  observed. 

A point  image  which  actually  covers  less  than  10  ^ of  the  active  area 
may  appear  to  cover  a substantial  amount,  10%  to  50%  or  even  all  of 
the  active  area,  if  it  is  bright  enough. 

Since  blooming  is  caused  by  charges  traveling  laterally  from 
diode  to  diode,  anti-blooming  provisions  can  be  incorporated,  which 
makes  it  very  probable  that  charges  will  be  absorbed  before  they  reach 
the  next  diode.  This,  however,  must  not  interfere  with  the  vertical 
travel  of  the  charges,  which  would  reduce  the  gain  of  the  target. 

The  most  effective  anti-blooming  targets  today  are  the  deep-etched 
devices  developed  by  Westinghouse  and  Texas  Instruments.  The  diodes  are 
separated  by  deep  furrows  which  prevent  the  depletion  layer  from 
extending  into  the  bulk  of  the  n-type  silicon  below  the  furrows.  A 
typical  curve  describing  the  blooming  characteristics  of  this  target  is 
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shown  in  Figure  III-l.  For  comparison,  a curve  with  a slope  of 
1/2  would  indicate  the  apparent  size  of  a point  source,  if  no  anti- 
blooming provisions  were  present  in  the  target. 

3.4  THE  READOUT  PROCESS 

3.4.1  Equivalent  Electrical  Circuit 

As  mentioned  in  the  preceding  section,  the  charge  carriers 
created  by  the  photoelectrons  in  the  bulk  n-type  silicon  partially 
discharge  the  capacitors  formed  by  the  diode  p-n  function.  In  other 
words,  the  reverse  bias  on  the  individual  diodes  is  reduced  from 
the  original  value  Vfi  to  Vc  (Figure  III-2) . The  electrons  in  the 
read-beam  "see"  different  potentials  and  by  landing  on  the  target 
restore  the  original  reverse  bias  on  each  diode.  The  number  of 
electrons  needed  to  recharge  the  individual  capacitors  is  equal 
to  the  number  of  carriers  created  by  the  photoelectrons;  these 
recharging  electrons  constitute  the  output  signal. 

To  analyze  the  readout  process,  it  is  necessary  to  draw  an 
equivalent  electrical  circuit  for  the  target  read  beam  interaction 

(Figure  III-2) . The  electrons  leave  the  thermal  cathode  and  are 
accelerated  by  the  voltage  V between  the  cathode  and  the  read  side 
of  the  target.  is  the  net  acceleration  voltage  for  the  electrons 

as  can  be  seen  in  Figure  II-6.  The  target  itself  can  be  considered 
to  consist  of  a capacitance  and  a diode  in  parallel  as  shown  in 
Figure  III-2.  A positive  bias  voltage  V„  (10  Volt)  is  applied 
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Figure  HI-I  TYPICAL  BLOOMING  CHARACTERISTICS  OF  EBSICON 
(DATA  FURNISHED  BY  S.  CAMPANA,  NAVAL.  AIR  DEV.  CENTR. ) 
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Figure  HE- 2 EQUIVALENT  CIRCUIT  FOR  READ  PROCESS 


to  the  n-type  silicon  bulk  (target  connector).  The  potential  dis- 
tributions in  the  target  as  they  change  during  read  and  exposure 
cycles  are  shown  in  Figure  III-3. 

The  situation  as  it  exists  before  any  voltages  are  applied  is 
shown  in  Figure  III-3a.  The  target  connector  is  at  cathode 
potential,  the  p-n  junction  causes  a small  negative  contact  potential 
Vcp  for  the  gold  cap  on  top  of  the  p-islands.  (Vcp^0.3  Volts,  the 
contact  potential  on  the  n*"  junction  is  also  approximately  0.3  V.) 

Figure  m-3b  shows  the  potential  distribution,  when  the  tareet  is 
connected  to  the  positive  bias  voltage,  V , but  before  the  first 

D 

scan  reaches  the  diode.  During  the  first  scan  the  p-islands  are 
placed  at  cathode  potential,  the  voltage  differential  across  the 
capacitor  becomes  now 

= V 

C B (3.4.1) 

This  voltage  differential  across  the  capacitor  occurs  mainly  in  the 
depletion  layer  of  the  p-n  junction,  while  the  rest  of  the  n and  p 
sections  have  basically  zero  potential  gradients.  The  depletion 
layer  increases  in  width,  depending  on  the  resultant  voltage 
(Figure  III-3c. 

During  illumination  of  the  target,  electron-hole  pairs  are 
created.  The  electrons  travel  towards  the  n+  layer  while  the  posi- 
tive holes  pass  through  the  diode,  discharging  the  capacitor  and 
thereby  reducing  the  voltage  Vp.  The  resulting  potentials  are  shown 
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Figure  HI-3  DISTRIBUTION  OF  POTENTIALS  ACROSS  STORAGE  TARGET 
(ONLY  UPPER  LIMITS  OF  CONDUCTION  BAND  ARE  SHOWN , CONTACT  POTEN- 
TIALS BETWEEN  METAL  CONTACTS  AND  SILICON  BULK  MATERIAL  ARE 
NEGLECTED) 
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in  Figure  III- 3d.  At  the  next  scan  the  read  beam  electrons  flow 
through  the  target  into  the  signal  connectors  until  the  voltage 
VT  has  become  zero  again  (Figure  III-3c) . 

3.4.2  The  Current  Voltage  Characteristics 

The  target  read  beam  interaction  is  quantitatively  determined 
by  the  so-called  current-voltage  characteristic  (I-V  curve) . This 
curve  gives  the  instantaneous  signal  current  I as  a function  of  the 
target  voltage  VT  for  a single  diode.  This  I-V  curve  and  its  particular 
shape  defines  in  principle  the  total  tube  performance,  the  signal  gen- 
eration in  the  tube,  and  the  tube  limitations.  This  I-V  curve  can  be 
determined  by  experiment  for  each  tube  under  operating  conditions. 

The  potential  VT  cannot  be  measured  directly  because  VT  is  the 
potential  difference  between  the  floating  electrode  of  the  storage 
target  and  the  thermal  cathode,  and  is  therefore  inaccessible  for  any 
experiment.  Indirectly,  however,  VT  can  be  obtained  by  saturating 
the  storage  target,  i.e.,  by  discharging  the  capacitors  in  the  p-n 
junction.  In  this  case  the  potential  Vc  across  the  capacitor  is 
reduced  to  the  known,  constant  contact  potentials  in  the  target. 

The  voltage  V^  now  becomes 


VT  = VB-  vT  (3. A. 2) 

This  important  1^  curve  can  now  be  determined  experimentally  by 
using  different  bias  voltages  Vg,  correcting  for  vT  and  measuring 
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the  corresponding  signal  current  I^.* **  A typical  I-V  curve  for  a 

silicon  diode  target  is  shown  in  Figure  III-4.  This  curve  has  all 

the  characteristics  of  a vacuum  diode  operating  with  a space 

charge  limited  regime.  Examination  of  Figure  III-4  shows  three 

distinct  areas:  the  region  of  exponential  increase  of  the  signal 

current  Ig  at  small  negative  potentials  VT;  the  space- charge  limited 

3/2 

region,  where  the  signal  current  increases  basically  with  , 

and  finally  a saturafcipn  region  where  the  signal  current  reaches 
the  value  of  the  saturation  currept  and  becomes  independent  of  VT. 

The  generally  assumed  fact,  that  the  target  simply  acts  as  an 
analyzer  of  the  energy  distribution  of  the  electrons  in  the  read 
beam, cannot  be  supported  by  the  experimental  evidence.  A detailed 
analysis  of  the  effects  determining  the  shape  of  the  I-V  curve  will 
be  the  subject  of  a separate  document  to  be  published  in  the  near  future. 

3.4.3  Approximation  of  the  I-V  Characteristic 

For  a mathematical  model  the  experimentally  obtained  I-V 
curve  can  be  approximated  very  closely  by  the  mathematical  form 

I = I (1  - e "b  VT)  <3*4-3> 

S ^SAT  u 6 ' 

where  b is  a constant. 

This  curve  gives  a good  representation  of  the  linear  range 

*Suggested  by  Dr.  A.  S.  Laponslcy  of  Westinghouse. 

**Plotted  from  data  furnished  by  R.  F.  Franseen,  Army  NVL. 
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Figure  HI-4  TYPICAL  I-V  CHARACTERISTIC  OF  EBSICON  TARGET 


and  the  transition  into  the  saturation  region.  Figure  III-5 
shows  both  curves,  normalized  for  equal  saturation  current.  Except 
for  the  region  around  the  origin,  both  curves  have  very  similar 
shape.  For  this  reason  a mathematical  model  should  give  very 
accurate  results,  except  for  very  low  lag  values,  which  will  occur 
only  after  many  reads  or  for  very  slow  scans.  But  these  effects  are 
considered  of  second  order  for  which  the  mathematical  model  was  not 
designed,  since  it  is  not  intended  to  operate  the  sensor  tubes  in 
slow  scan  or  with  many  read  scans  in  a real  system  application. 

3.4.4  The  Signal  Current 


As  observed  from  a fixed  point  on  the  target,  the  passing  of 
the  read-beam  appears  to  be  equivalent  to  the  arrival  of  a time- 
varying  current.  If  the  beam  scan  velocity  is  constant,  the  time 
dependence  of  the  current  follows  a curve  similar  to  the  character- 
istic which  describes  the  spatial  current  distribution  in  the  read- 
beam.  If  this  distribution  is  given  by  a function  g(t)  and  the 
saturation  current  I is  replaced  by  the  total  read  beam  current 
1^,  Equation  (3.4.3)  can  be  rewritten 


1 


s 


(t) 


g(t)  (1 


(3.4.4) 


Since  is  determined  by  the  charge  Q stored  on  the  target  and  the 
current  I^(t)  is  defined  by 


I (t)  = (3.4.5) 

S v ' dt 

a differential  equation  for  the  signal  current  I^(t)  can  be  derived. 


42 


SIGNAL  CURRENT 
(ARBITRARY  UNITS) 


1 i 


lo 


0 

Vt=Vb-^t  TARGET  VOLTAGE 

( ARIBITRARY  UNITS) 


Figure  HI-5  APPROXIMATION  OF  EXPERIMENTALLY  OBTAINED 
I-V  CHARACTERISTIC  BY  FUNCTION  Is  = I0  ( I -e  - bVT) 
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The  detailed  derivation  and  solution  of  this  differential  equation 
for  a two-dimensional  charge  field  and  a two-dimensional  read  beam 
moving  with  constant  read  velocity,  v,  are  presented  in  Section  4.1 
of  this  report. 

The  solution  for  the  time  dependence  of  the  signal  current  for 
a given  charge  field  is  involved  and  requires  the  use  of  computers. 

3.4.5  The  Lag 

The  charge  stored  on  the  target  is  not  read  out  fully  during 
one  scan  process;  this  can  have  two  reasons*  both  of  which  limit  the 
signal  current. 

a)  Read  beam  limitation.  In  this  case  the  scan  beam  cannot 
furnish  enough  electrons  during  its  "dwell  time"  to  recharge  the 
target  and  a certain  percentage  of  the  original  signal  remains  on 
the  target.  If  there  is  true  beam  limitation  the  signal  current 
remains  constant,  independent  of  the  original  illumination;  the 
remaining  charge  is  merely  the  difference  between  the  charge  created 
by  the  illumination  and  the  total  charge  delivered  by  the  full  read 
beam  during  the  "dwell  time". 

b)  Lag . As  the  I-V  curve  shows,  the  signal  currents  vanish 
for  small  voltages  V^.  As  soon  as  the  target  potential  reaches 
cathode  potential,  less  and  less  electrons  are  available  due  to  the 
space-charge  effects  and  cannot  land  on  the  target.  This  effect  can 
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be  considered  a type  of  read-beam  limitation,  but  here  the  available 
charge  is  integrated  not  over  the  whole  read  beam  delivered  during 
the  "dwell  time",  but  only  over  the  number  of  electrons  in  the  read 
beam  having  an  energy  higher  than  a certain  cutoff  energy. 

Both  read- beam  limitation  and  lag  are  strongly  dependent  on 
the  scan  velocity,  since  this  is  the  parameter  which  defines  the 
"dwell  time".  Their  dependence  on  the  original  illumination  level, 

however,  shows  a different  behavior:  In  case  of  read-beam  limitation 

a constant  amount  of  charge  QgAT  is  replenished  during  a scan, 

independent  of  signal  level,  while  in  the  case  of  lag  a constant 

amount  of  charge  Q is  left  on  the  storage  target  during  each  scan, 

£ 

independent  of  signal  level.  Since  lag  values  are  usually  given  as 
a fraction  e of  the  actually  read  charge  to  the  amount  of  charge 
stored  on  the  target  Qq,  one  obtains  for  read-beam  limitation: 

Qo  ” QSAT  (3.4.6) 


and  for  lag 


e = . (3.4.7) 

L Qo 
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SECTION  IV 


THE  MATHEJ1AT I CAL  TUBE  MOPEL 

The  objective  of  this  section  is  to  establish  a representative 
model, and  to  do  so  the  qualitative  results  of  Section  3. A must  be 
carried  to  their  mathematical  conclusions  in  the  form  of  equations 
yielding  quantitative  values  of  readout  current  and  instantaneous  charge. 
Since  these  equations  are  functions  of  initial  target  charge,  the  specific 
probability  density  functions  of  the  processes  described  in  Sections 
3.1  through  3.3  which  describe  the  charging  of  the  target  must  be 
derived  and  justified.  This  is  accomplished  in  Section  A.l  below. 

In  Section  A. 2 these  results  are  brought  together  in  the  forms 
of  a computer  simulation  of  the  tube  operation  and  as  proprams  used 
to  evaluate  the  statistical  performance  of  the  tube. 

A.l  DESCRIPTION  AND  DEFIVATION  OF  EQUATIONS 

A. 1.1  The  peadout  Process 

In  Section  3. A,  the  relationship  between  diode  current  and 
voltage  and  hence  also  for  the  distribution  of  target  current  density, 

J^,  beam  current  density,  J^,  and  target  voltage  was  shown  to  be 
closely  approximated  by: 

-bV 

Jt  = Cl  -e  E)Jb  (4.1.1) 

2 

With  knowledge  of  target  capacitance  (farads/cm  ),  (A. 1.1)  can  be 
rewritten  in  terms  of  target  charge  density,  q^_,  as 


A6 


Jt  = (1  - e 


9 


(4.1.2) 


-gq 


t 


> Jb 


where , 


(4.1.3) 


The  readout  beam  is  assumed  to  be  a moving  axially- symmetric 
Gaussian-shaped  current  distribution  vrith  a total  current  equal  to 
as  follows: 

[x-x0(O]2 

T , N 0 

Jb  (x,y,t)  = — - e 

TTCC 


+ ,,  , 

— O (^.1. 


where  (x,  y)  is  a specific  point  on  the  target  and  Tx^(t),  y^t)] 
is  the  present  location  of  the  center  of  the  beam.  The  forms  of 
Xg(t)  and  YqCO  are  such  as  to  achieve  the  raster  scan-pattern 
described  in  Section  2.3. 

The  distribution  of  the  target-charging  current  density  is  then, 


dq 


t 


dt 


J 

t 


(4.1.5) 
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There fore, 


dqt(x,y,t) 

dt 


-6q  (x,y,t)  -l 

6 J Jbkc»y»t^ 


(4.1.6) 


suppressing  the  notation  x and  y and  rearranging  terms,  (4 . 6 . 5) 
becomes 


t 

- J Jb  (r)  d r - 

O 


Aq  (x,y,t) 


(4.1.7) 


The  negative  of  the  integral  of  the  beam- current  distribution,  Aq, 
will  be  retained  inexplicitly , for  now,  because  it  will  take  on 
different  forms  depending  upon  how  it  is  used  in  f^nal  implementations. 

Equation  (4.1.7)  can  be  integrated  to  solve  for  charge  density 
and  then  for  the  target  current  distribution  yielding  the  following: 


qt(t)  = ? In 


BAq(t)  r 6qt(0) 
e e 


-1] 


+ 1 


(4.1.8) 
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Then, 


J 

t 


(t)  = - 


dqt(t) 

dt 


8Aq(t)  r 8q  (0)  -| 

e Li 

8Aq(t)  r 8q  (0)  -| 

e L6  " l|  + 1 


(4.1.0) 


Here,  J,  (t)  also  is  not  explicitly  given  since  its  specific  form 
b 

depends  upon  the  implementation  of  these  results. 

Equations  (4.1.8)  and  (4.1.9)  are  general  closed-form  results 
for  any  target  with  the  J-q  characteristics  of  (4.1.2).  However, 
for  the  diode-array  target  of  an  EBSTCON  tube  these  equations  must 
be  integrated  over  the  area  of  a diode  so  that  the  individual 
current  contributions  from  each  diode  to  the  total  target  current 
can  be  determined. 

Since  each  diode  has  a gold  cap,  as  illustrated  in  Section  1.3, 
the  charge  and  its  derivative  are  constant  over  the  area  of  each 
diode.  It  will  be  assumed  that  the  read  beam  is  also  constant 
over  a diode  since  its  area  is  small  compared  to  the  read  beam  area. 
Hence,  for  each  diode,  (4.1.8)  and  (4.1.9)  become 


Q (x1,yj,t)  = J In 


8AQ(x_^,y.  ,t)  p BQ(x^,y.,o) 
J J e J 


+ 1 


(4.1.10) 
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and 


it(xj  .y^t) 


eBAQ(x±ty  ,t) 
e3AQ(xi,y^ ,t) 


feBQ(x^  ,Yj,0) 

[eBnCxi,yi>o) 


-t]  yy^t) 

-1]  +1 


(4.1.11) 


where 


AQ(t)  = Aq(t)  AA. 


(4.1.12) 


It(t)  = Jfc(t)AA 


(4.1.13) 


and 


0(t)  = qt(t)AA 


(4.1.14) 


0(0)  = qt(0)AA 


(4.1.15) 


in  all  of  which  AA  is  the  area  of  a diode  and  in  calculating  C in 
(4.1.3)  c must  also  be  multiplied  by  AA.  The  read-beam  current  is 
now 


Ib(xi,yj,t) 


T0AA 


[xi-xQ(t)]2  + [y.j-y0(t)]2 

2 

a 


(4.1.16) 
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where  (x 


is  the  center  of  the  i,j 


th 


diode . 


*t  ^Xi,yj  iS 


i 


the  independent  contribution  to  the  total  target  current  by  the  i, 


determine  the  target  current  at  tine  t. 

Equations  (4.1.10)  and  4.1.11)  provide  a deterministic 
functional  relationship  for  charge  and  current  contribution  for  each 


4.1.2  Statistics  of  Photoelectrons  Emitted  by  the  ^hotocathode 

If  statistical  analyses  are  to  be  performed,  the  probability 
density  function  (PDF)  of  the  initial  charge  must  be  known.  There- 
fore, in  this  and  the  following  two  subsections,  the  processes 
described  in  Sections  3.1  thru  3.3  must  be  analyzed  to  derive  and 
verify  the  specific  PDF's  of  those  processes  leading  to  the  final  PDF 
required;  namely  that  of  0(x^,y^,0). 

It  is  shown  in  [1,  P.  496]  that  for  the  case  where  the  observation 
time,  T,  is  large  compared  to  the  coherence  time,  x,  the  probability 
of  the  number,  N,  of  photoelectrons  emitted  to  each  diode  in  the 
time  interval  T is: 


diode  and  therefore  they  must  be  summed  over  all  i and  j to 
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where 


6 


(1  + 2 N 1)1/2 


(4.1.17) 


and 


so(f)= 


si<T>  = 1 


(4.1.18) 


but  for  n greater  than  one,  from  [2]  and  [4,  P.  467], 


S ( 


n 


n-1 

i + E 

k=i 


(n-l-fk)! 
k'.  (n-l-k)  !2k 


(4.1.10) 


In  the  case  of  the  EBSICON  (Section  3.1),  the  observation  tine, 

T (integration  time)  is  much  longer  than  the  coherence  time,  x, 

-13 

of  the  light  illuminating  the  photocathode  (x/T  « 1.5  x 10  ). 

In  this  case,  6 approaches  unity  so  that 


Prob . 


{N=n}  = 

n. 


(g-D, 

t/T  J 


(4.1.20) 


However,  the  exponent  approaches  the  indeterminate  form  0/0.  But, 
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using  LfHospitalfs  rule,  it  is  easily  shown  that  this  form  is  resolved 

and  the  limit  becomes  -N.  Also  except  for  values  of  n in  the  order 

of  10^,  corresponding  to  extremely  great  illumination  (ordinarily 

3 

n ranges  from  50  to  10  ) , it  is  seen  from  (4.1.19)  that  Sn(T/x) 
is  very  nearly  unity.  Therefore,  (4.1.20)  approaches  the  Poisson 
probability , 


Prob . 


N 

{N=n}  = 

n. 


(4.1.21) 


4.1.3  Gain  and  Loss  Mechanisms 

As  described  in  Section  3.2,  gain  is  achieved  in  the  tube  by 
virtue  of  accelerating  the  photoelectrons  emitted  from  the  photo- 
cathode in  the  image  section,  whereupon  they  bombard  the  target  and 
create  hole-electron  pairs  in  the  n-type  silicon  layer.  The  number 
of  hole-electron  pairs  created  by  each  photoelectron  determines  the 
gain. 

Since  all  photoelectrons  arrive  with  essentially  the  same 
energy, and  since  in  giving  up  energy,  while  being  absorbed,  it  is 
given  up  in  equal  discrete  amounts  (each  amount  corresponding  to  the 
energy  required  to  create  a hole-electron  pair,  3.4eV),  the  gain 
process  is  deterministic  for  all  practical  purposes  for  a constant 
rate  of  arrival  of  photoelectrons.  Further  since  the  cross  section 
of  the  target  is  high,  the  holes  created  proceed  under  a slight 
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field  gradient  toward  the  p-islands  very  evenly  spaced  in  time  while 
the  electrons  remain  in  the  n- region. 

Random  loss  occurs  through  recombination  of  some  of  the  hole- 
electron  pairs  in  the  n-silicon  layer.  If  a given  interval  in  time 
is  considered  as  the  holes  travel  in  an  evenly- spaced  stream,  each 
hole  is  either  recombined  or  not  with  some  probability,  p.  There- 
fore, the  number  of  holes  remaining  in  a given  interval  of  time  as 
they  arrive  at  a p-island  is  governed  by  a Bernoulli  distribution. 
This  is  because  the  success  or  failure  of  the  arrival  of  each  hole 
in  the  interval  was  determined  by  an  independent  Bernoulli  trial. 

It  is  stated  in  Section  3.2,  however,  that  these  losses  are 
low.  This  statement  is  equivalent  to  the  statement  that  the  proba- 
bility, p,  is  approximately  zero  and  that  the  standard  deviation 
of  the  Bernoulli  distribution  is 

a = N/np(l-p)  = 0 (4.1.21) 

thus  yielding  a process  which  is  essentially  deterministic. 

The  above  discussion  was  based  upon  a constant  rate  of  arrival 
of  photoelectrons  when,  in  fact,  the  rate  of  arrival  was  shown  in 
Section  4.1,2  to  be  Poisson  distributed.  Therefore,  in  the  case  of 
the  EBSICON  tube  the  Poisson  statistics  of  photoelectron  count  are 
passed  over  to  stored  charge  by  the  tube  gain  factor,  G, 
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The  resultant  distribution,  which  can  be  called  a ,Jmodif iedn 
Poisson  distribution  [3,  P.  33]  is  obtained  by  multiplying  a Poisson 
distributed  random  variable,  x,  by  a constant,  G,  to  yield  a new 
random  variable,  y. 

y = Gx  (4.1. 22) 

If  the  mean  of  the  distribution  of  x is  A,  its  PDF  is  given 
by: 

°° 

fx(*)  = £ \r-  <5(x-k)  (4.1.23) 

k=0 


and  the  PDF  of  y is 


00  e-*xk 

f (y)  = L - 7-,  6(y-Gk)  (4.1.24) 

y k=0  K* 

with 

V = GX 

y (4.1.25) 

and  o2  = C2X  (4.1.26) 

y 

4.1.4  Target  Saturation 

As  described  in  Section  3.3,  there  is  local  diode  saturation 
with  high  illumination  which  in  turn  leads  to  blooming.  If  the 
random  variable,  y,  in  (4.1.24)  is  the  available  or  total  possible 
charge  per  diode  with  no  saturation  effects,  then  the  statistics  of 
the  actual  charge,  Q^,  are  changed  from  those  of  y by  saturation. 
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(4.1.27) 


The  functional  relationship  chosen  between  y and  to 
describe  local  saturation  is 

Q.  = ± (1-e'YY) 

c Y 

This  function  was  chosen  for  four  reasons.  The  first  two  are  the 
facts  that  in  the  linear  region  where  saturation  is  not  significant, 
i.e.,  where 


-x  'v  , (4.1.28) 

e « 1-x 


the  mean  and  variance  of  Qt  can  be  shown  to  approach  the  mean  and 
variance  of  y,  given  in  (4.1.25)  and  (4.1.26)  respectively. 

The  other  two  reasons  are  that  when  saturation  occurs  (i.e., 
when  X-^o°)  the  mean  of  Q approaches  a constant,  1/y,  and  the  variance 
of  Q approaches  zero.  These  two  effects  are  the  most  significant 
features  of  local  saturation  and  are  clearly  seen  in  the  following 
expressions  (4.1.29)  and  (4.1.30)  for  the  mean  and  variance  of  Q^_ 
respectively. 


A n-,-x<1-e‘YG> 


\ ’ Y [1-6 


(4.1.29) 
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-X[l-(e  YG)21  -2X(l-e"YG) 

° “6 


(4.1.30) 


The  transformation  of  (4.1.27)  yields  the  following  probability 
function  (PDF)  for  Qt 


-A  k 
e A 


«[Qt-  i (l-e'^Gk)] 
Qt  c k=0*C-  fc  Y JJ 


(4.1. 31) 


This  PDF  for  the  initial  diode  charge  is,  therefore,  the  one  to  be 
used  in  computer  implementations  which  evaluate  the  statistical 
performance  of  the  tube.  Blooming  is  accommodated  in  the  imple- 
mentations by  the  proper  initialization  of  the  charge  of  the  diode 
array,  as  will  become  evident  in  Section  4.2. 

4.2  IMPLEMENTATION  IN  TERMS  OF  COMPUTER  PROGPAMS 

Equations  (4.1.10),  (4.1.11)  and  (4.1.31)  are  the  fundamental 
building  blocks  in  implementing  the  mathematical  model  as  derived 
-in  the  preceding  section.  They  can  be  used  in  two  ways.  One  way 
is  to  simulate  the  video  signal  or  the  sampled  video  output  of 
the  tube.  The  other  way  is  to  determine  the  ensemble  variance  and 
mean  at  any  point  on  the  "current  map"  of  the  output  and  the  correla- 
tion between  two  separated  points  on  the  current  map.  From  these 

results  P and  P can  be  determined  for  various  system  conditions 
D r A 

for  a given  threshold.  Computer  programs  were  written  to  accomplish 
both  implementations . 
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A. 2.1  Simulation  Program 


In  the  simulation  program  the  mean  number  of  photoelectrons 
arriving  from  the  photocathode  during  the  integration  time,  A^  , 
is  determined  for  each  i,^*1  diode.  This  number  is  dependent  upon 
the  optics  assumed  for  the  system,  the  background  level  assumed  and 
the  strength  of  the  image,  if  any,  illuminating  the  i,:)^  diode. 

The  Poisson  parameter  is  used  in  the  following  equation  to  approximate 
a Poisson  random  variable  with  a Gaussian  envelope,  but  taking  on 
integer  values: 


n^  = Integer 


/-2A_. .In 


ij±nuL 


cos(2iru_) 


+ A , 


ij 


.5 


(4.2.1) 


where  u^  and  u^  are  two  independent  pseudo- random  variables  which  are 
uniformly  distributed  over  the  interval  ranging  from  zero  to  one. 

The  samples  of  u^  and  u^  are  generated  by  means  of  a power  residue 
method  for  each  diode. 


From  (4.1.27)  the  charge  for  each  diode  is  computed  by 


q = — (i  _ e-GYnlj)  (4.2.2) 

ij  Y ' 

where  — is  the  maximum  or  saturation  value  for  target  charge  per 
diode  and  G is  the  tube  gain,  multiplied  by  the  electronic  charge. 

In  this  application, ^Q(x^ ,y^ (t)  of  (4.1.12)  is  appropriately 
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I AAe 
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/Z£ If 
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/X  -X  -vt  \ 

erf(  — ) (4.7.3) 
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-erf 


where  L is  the  distance  between  scan  lines  and  n is  the  line  number 
of  the  present  scan.  The  scanning  bean  is  as  defined  from  (4.1.1b) 
by  letting  x^(t)  = -xo~vt  w^ere  v the  scan  velocity  and  be  lettling 
y^(t)  = Q^.  From  (4.1.8),  0(x_^,y_.,t)  can  be  calculated  for  any 
tine,  t,  along  the  present  scan  line  for  each  diode. 


Q(x  ,y  ,t)  = 1 In' 
1 J ^ 

eBAQ(xi5y^ ,t) 

e3Q(xi,yj,0)_1 

J 

+ 1 

(4. 2. A) 


For  the  first  scan  line  Q(x_^,  y^,0)  is  given  by  (4.2.2)  and  for 
each  scan  line,  t varies  from  zero  to  -2^/v  where  Xq  is  a large 
negative  distance  relative  to  the  diode  array  size. 


In  the  particular  sinulation  an  array  of  21x21  diodes  was  used 
with  a spacing  between  diode  centers  of  12.5ym  or  .00125  cn. 


The  origin  of  the  coordinate  system  of  x^  ,y_.  of  the  array  was  taken 

to  be  the  center  diode  (i.e.,  • Therefore,  x and  y ranged  from  -.0125 
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to  +.0125  cm  and  x was  set  at  -.025  cm. 

0 

After  calculating  the  new  charge  for  each  t on  the  first  scan 
line,  the  individual  diode  current  contribution  is  calculated  as 


Mxi,yj,t:^  = 


e-6Q(xi#yj ,t) 


Vw0 


(4.7.5) 


Then  for  each  t the  total  target  current  is  calculated  as 


21  21 


:(t)  = E E ir(x.>y. >t) 

i-1  .1=1  J 


(4.2 .ft) 


At  the  end  of  each  scan  (i.e.,  when  t = -2xq/v)  , the  old  0(x^,y^,0) 


is  updated  as 


Q(x.,y.  ,0)  = Q(x  ,y 2-  ) 

1 J 1 J V 


(4.2.7) 


This  value  is  then  used  in  (4.2.4)  to  calculate  the  charge,  current 
and  final  charge  for  the  next  scan.  The  above  procedure  is  repeated 
until  the  last  scan  is  completed. 

If  the  ‘'increments  in  time"  are  made  fine  enough,  a three- 
dimensional  plot  of  the  simulated  video  output  such  as  shown  in 
Figure iv-1  can  be  made.  Note  that  the  time  axis  (horizontal  axis) 
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Figure  IV-1  SIMULATED  VIDEO  OUTPUT 


is  labeled  in  centimeters  and  hence  the  use  of  the  term  ’current 
map'’.  Shown  in  Figure  IV-1  is  the  current  map  for  background  with  two 
point  images  separated  by  four  diodes.  On  the  other  hand,  one  may 
want  to  make  the  time  increments  correspond  to  sampling  intervals 
to  be  used  for  video  sampling.  Then  the  simulation  cou]d  be  used 
to  print  out  just  the  values  of  the  sampling  points. 

4.2.2  Statistical  Evaluation  Programs 

The  main  idea  behind  the  programs  which  calculate  the  ensemble 
statistical  parameters  is  to  first  integrate  the  read  beam  to  a point 
(in  the  specific  case,  to  the  center)  in  the  current  map  and  to  write 
the  equation  for  the  beam  current  evaluated  at  that  point  also.  A 
13  x 13  diode  array  was  chosen  in  this  case  to  ensure  that  edge 
effects  would  not  be  a factor  in  computing  the  output.  In  this 
case,  AQ  and  I would  be  calculated  as  follows: 

(4.2.3) 
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(A. 2.0) 


VW  = ZT  e 


x.2  + y.2 
i 3 


a*1 


TTOt 


The  next  thing  to  recognize  from  (4.1.11)  is  that  the  tareet 
current  contribution  from  each  diode  is  a function  of  its  initial 
charge  written  in  a closedform  expression: 


0AQ(x  ,y  )T  8Q(x.,y  0) 

= e J Le  1 ] -1 

8AQ(x  y )[  8Q(x  y 0)  f 
e JLe  J -1J  + 1 


(A. 2.10) 


Then  with  the  probability  density  function  ("PDF)  given  in  (A  .1.31) 
for  each  diode, 


vv 


- E 

k=0 


k! 


Tq. . - ~(1  - 

I ij  Y 


-kyG. . 
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J 


(4.2,11) 


the  mean,y  (x  , y . ),  and  the  mean  square,  (x^jy^),  of  the  current 
It  1 1 t 1 

contribution  for  each  diode  can  be  closely  approximated  as 


(x 


i'V 


= Z 

k=0 


k 

e Zhl! 
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(4.2.12) 
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(4.2.13) 


Since  the  weights  in  the  summations  are  the  Poisson  weights  of  the 
original  Poisson  density  of  the  number  of  photoelectrons,  k^  was 
chosen  to  be  equal  to  an  integer  rounded  to  the  mean  plus  three 
times  the  standard  deviation  of  that  density, 


= Integer 


[*u + 


vA 


ij 


+ 


(4.2.14) 


From  (4.2.12)  and  (4.2.13)  the  variance  for  each  diode  current  can 
be  calculated, 


(Xi ,yj  > = ml  ^xi 


•yi>  "“i 


(4.2.  !■=;) 


Each  contribution  is  statistically  independent  of  all  the  others  and 
hence  the  mean  and  the  variance  of  the  total  target  current  can  be 
obtained  by  summing  the  means  and  variances  of  all  the  contributions, 


13  13 

yT  ■ L L uT  (x. ,y  ) 
1 i=i  j-i  t 1 J 


(4.2.16) 
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and 


°T2  - L L a 2 (x  y ) 
i=l  j=l  1 J 


(4.2. 17) 


To  calculate  the  output  after  the  next  write-read  cvcle  the 
mean  and  variance  of  the  charge  remaining  on  each  diode  after  a 
complete  read-out  must  be  computed.  Since  the  charge  from  the 
next  illumination  is  independent  of  the  remaining  charge,  the  two 
means  and  variances  are  summed.  Then  from  the  new  variance  and 
mean,  the  corresponding  A and  G_^  must  be  calculated  for  use  in 
(4.2.11).  The  integrated  beam  current  AO(x_^,y\)  must  correspond  to 
the  final  time  after  a complete  read-out. 


and  again  recognizing  from  (4.1.10)  that  the  final  charge  is  a 
function  of  the  initial  charge 


(4.2.18) 


2 


AQ  (x.,yj)  = 


Qf  xi>yj»  Q (x^y.,0) 


(4.2.1?) 


= « In  e 
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the  mean,  mean  square  and  hence  the  variance  can  be  found  in  the  same 
manner  as  these  parameters  were  found  for  each  diode  current  (A. 2. 12), 
(4.2.13)  and  (4.2.15). 


For  the  nev;  illumination,  the  mean  charge  and  the  variance  can 
be  found  from  (4.1.29)  and  (4.1.30)  respectively, 


VW  - 7 


-A  (1  - e-Y°) 
1 - e 


(4.2.20) 


and 


a 2 
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(x. ,y . ) = - 


i j 


r-\  a 

e J 


(1 


(4.2.21) 


To  (4.2.20)  is  added  the  mean  of  the  remaining  charge  on  the  i,  diode 

and  to  (4.2.21)  the  variance  of  the  remaining  charge  is  added.  Let 

the  resultant  mean  and  variance  be  designated  by  pQ(x^,y  ) and 
2 

o p(x^,y^).  Then  by  inverting  (4.2.20)  and  (4.2.21)  simultaneously 

the  appropriate  A.,  and  G..  to  be  used  in  (4.2.11)  can  be  calculated 

ij  1.1 

by 


Jln[l-YpQ(Xi,y .)] 

In 

i-Yy0(xi,y1) 
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4.2.22 
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and 


C = - - In 

ij  y 


ln|y2aQ  (x1»yi)+[1  ~ YyQ  (Xi  »yi  )]  2 1 
ln[l  - YPQ(x.,y.)J 


-1 


(4.2.23) 


The  above  steps,  equations  (4.2.8)  through  (4.2.23),  are  repeated 
for  each  subsequent  read-out  and  illumination  cycle. 

Also  included  in  the  programs  is  the  ability  to  calculate,  in 
a similar  manner,  the  means  at  two  different  points  in  the  current 
map  and  the  expectation  of  the  product  of  the  currents  at  two 
different  points.  Of  course,  AQ,  and  hence  I will  be  different 
at  each  point  but  obtaining  the  forms  is  straightforward  and  will 
not  be  explained  here.  However,  it  must  be  stated  that  in  order  to 
use  the  PDF  of  (4.2.11)  and  procedures  such  as  those  of  (4.2.12) 
and  (4.2.13)  the  illumination  must  be  stationary  background  noise 
for  each  cycle.  That  is  to  say,  A^  must  be  the  same  for  the 
charge  at  the  two  points.  This  is  not  the  case  when  only  one  point 
is  used.  In  fact,  when  one  point  is  analyzed,  the  procedure  allows 
a scene  to  be  illuminated  and  for  the  scene  to  change  from  cycle  to 
cycle*  since  A_^  can  be  chosen  by  the  user  to  be  different  for  each 
i and  each  j. 
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Finally,  from  a separate  analysis,  the  PDF  for  an  individual 
diode  current  contribution  was  determined  and  it  was  found  to  closely 
match  a Gaussian  PDF.  Therefore,  the  PDF  of  the  sum  of  all  I^x^y^) 
will  be  even  more  Gaussian.  Hence,  by  knowing  o and  y , the  PDF  of 
the  target  current  is  known.  By  analyzing  one  point  first  without 
an  image  present  but  with  only  background  noise  and  then  with  an  image 
in  the  presence  of  background,  the  two  PDF’s  can  be  superimposed, 
thereby  allowing  P^  and  P^  to  be  determined  for  a given  threshold. 
Further,  by  analyzing  two  points  while  increasing  the  distance  between 
them  on  each  analysis,  the  autocovariance  function  of  the  background 
noise  can  be  found.  Then  by  taking  the  Fourier  transform,  the  noise 
power  spectrum  is  known  and  from  it  the  noise  bandwidth  can  be  determined. 
4.2.3  Computer  Listings 

The  listings  of  the  programs  described  in  Section  4.2  are 
provided  in  the  appendices  in  the  following  order:  The  simulation 

programs  (Appendix  A)  and  the  two  statistical  evaluation  programs 
(Appendices  B and  C) . The  last  two  listings  are  essentially  the 
same  except  that  the  first  allows  the  computation  of  the  normalized 
auto  covariance  function  of  background  noise  whereas  the  second  does 
not.  The  second  program  was  specifically  written  for  the  inclusion 
of  images  as  well  as  background  noise. 

The  specific  parameter  values  listed  and  the  arrangement  of 
DO-loops  are  changed  in  all  the  programs  to  accommodate  various 
operating  conditions.  As  listed,  only  a specific  tube  and  operating 


mode  is  represented. 
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SECTION  V 


RESULTS  AND  COMPARISONS  WITH  EXPERIMENTAL  DATA 

The  simulation  program  described  in  subsection  A. 2.1  and  the 
statistical  evaluation  program  described  in  subsection  A. 2. 2 were 
run,  and  when  available,  experimental  data  were  used  to  examine 
the  trend  results  of  the  model.  In  some  cases,  specific  tube  and 
operating  parameters  were  known  and  were  inserted  into  the  model  so 
that  direct  instead  of  trend  comparisons  could  be  made. 

Of  interest  were  flat-field  response  (i.e.,  only  background 
noise  was  present),  point-source  response,  response  to  images  in 
the  presence  of  background  noise,  and  resolution.  Within  these 
broad  interests  were  some  critical  interest  items  of  tube  response. 
These  included  saturation  effects,  image  splitting,  point  spread  and 
modulation  transfer  functions,  signal  to  noise  ratio,  probability  of 
target  detection,  and  probability  of  false  alarms  caused  by  back- 
ground noise  for  given  threshold  settings. 

5 . 1 Flat  Field  Response 

The  signal  current  as  a function  of  uniform  illumination  of  the 
photocathode  is  defined  as  the  flat  field  response.  In  Figure  V-I, 
typical  response  curves  are  shown  for  an  EBSICON  tube.  These  curves 
are  based  on  measurements  made  by  the  ITEK  Corporation  on  an  RCA 
A80A/SUR  sensor.  The  top  curve  shows  the  flat  field  response  of 
the  mean  signal  current  while  the  lower  curve  shows  the  behavior  of 
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the  RMS  noise  values  associated  with  the  mean  signal  currents.  The 


flat  field  response  clearly  has  two  distinct  sections:  The  straight 

line  inclined  at  45°  representing  the  linear  part  of  the  sensor 
response  ranging  over  about  two  decades,  and  the  horizontal  part  of 
the  curve  indicates  saturation  of  the  sensor.  The  curve  represent- 
ing the  RMS  values  also  shows  several  interesting  features.  For 

both  very  low  and  very  high  input  levels  this  curve  approaches  a 

-9 

constant  in  the  vicinity  of  3 x 10  amp.,  the  preamplifier  noise. 

If  this  preamplifier  noise  is  subtracted*  from  the  measured  RMS 
values,  a straight  line  is  obtained  with  a slope  of  one  half.  This 
curve  represents  the  photo-electron  noise  limit.  At  the  saturation 
point  the  measured  RMS  value  drops  rapidly  to  the  amplifier  noise. 
This  corresponds  exactly  to  the  expected  behavior. 

These  experimental  results  were  used  for  comparison  with  the 
output  of  the  computer  model  (Figure  V-2) . The  curves  shown  in 
Figure  V-2  include  those  of  Figure  V-l,  replotted  with  the  pre- 
amplifier noise  subtracted.  Four  curves  are  shown;  the  solid  lines 
represent  the  computed  values  and  the  dotted  lines  are  the  experi- 
mental results.  Both  sets  of  results  were  obtained  independently 
and  several  operational  parameters  were  different,  and  for  this 
reason  their  influence  must  be  considered  in  order  to  make  a 

^Subtraction  here  is  understood  to  apply  to  the  variances  of  both 
the  measured  values  and  the  preamplifier  noise 


Measured  -a  Preamplifier 
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Figure  3Z-2  EXPOSURE  ON  PHOTOCATHODE 
(PHOTOCATHODE  NORMALIZED  TO  BE  EQUAL  TO  TARGET  SIZE) 
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quantitative  comparison.  But  some  basic  facts  can  be  obtained 


immediately : 

o The  computed  flat  field  response  curve  has  the  same  shape 
as  the  experimental  one;  a linear  range  (unity  slope)  and  a 
sharp  knee,  where  it  turns  into  the  saturation  region 
(horizontal  line) . 

o The  curve  representing  the  RMS  values  again  has  the  same  shape 
as  the  experimental  curve;  it  follows  a line  with  a slope  of 
one  half  and  at  the  saturation  point  it  drops  to  zero  very 
rapidly. 

In  order  to  compare  the  curves  quantitatively  the  following  facts 
have  to  be  considered: 

o The  areal  scan  speed  in  the  case  of  the  computer  simulation 

was  9 x 10^  mm^/sec,  but  in  case  of  the  experiments  it  was 

3 2 

9.4  x 10  mm  /sec.  Since  the  signal  current  in  case  of  flat 
field  illumination  is  directly  proportional  to  the  areal 
scan  velocity,  it  would  be  expected  that  the  computer  simu- 
lation values  would  be  higher  than  the  experimental  results 
by  a factor  10 . If  the  experimental  curve  in  Figure  V-2  is 
shifted  by  this  amount,  the  linear  sections  of  both  curves 
coincide  exactly. 

o In  the  computer  simulation,  read  beam  saturation  took  place, 
while  in  the  experiments,  target  saturation  occurred.  The 
read  beam  limitation  for  the  experimental  curve  should  be 


73 


approximately  6 to  7 x 10  ^ amp.  If  the  read  speed  in  the 
experiments  had  been  increased  by  a factor  10,  beam  limitation 
would  have  occured  at  this  value  and  both  the  experimental 
and  computed  response  curves  would  coincide  entirely, 
o The  RMS  noise  current  is  proportional  to  the  square  root  of 
the  areal  scan  velocity;  in  this  case  the  computed  values 
would  be  expected  to  be  higher  than  the  experimental  results 
by  a factor  of  approximately  3.3.  Figure  V-2  shows  that 
the  factor  is  approximately  2.7.  This  discrepancy  can  be 
explained  by  correcting  for  the  smaller  read  beam  diameter 
used  in  the  experiments,  which  results  in  a smaller  integration 
area  and  therefore  larger  noise  values.  If  the  experiments 
had  been  run  under  beam  limiting  conditions  with  10  times 
the  scan  velocity,  both  RMS  curves  would  coincide  also. 

In  other  words,  by  properly  accounting  for  the  different  operational 
parameters,  the  curves  can  be  normalized  and  the  resulting  pairs  of 
curves  in  Figure  V-2  could  be  made  to  coincide  in  slope,  absolute 
values  and  in  saturation  values.  The  comparative  results  are  encour- 
aging. Additional  simulations  with  matched  parameters  are  planned 
and  should  furnish  further  verification  of  the  mathematical  model. 

5.2  POINT  SOURCE  RESPONSE 

Shown  in  Figure  V-3  is  the  basic  point  response  characteristics 
of  the  model  obtained  from  the  simulation  program  described  in 
section  4.2.1.  The  response  is  presented  in  the  form  of  a signal 
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Figure  V-3  POINT  RESPONSE 


current  "map"  which  geometrically  overlays  the  diode  array  (not 
shown) . The  read  process  goes  from  left  to  right  and  from  bottom 
to  top.  Several  effects  are  to  be  observed;  the  point  readout  char- 
acteristic is  asymmetric  in  both  directions  and  there  is  image 
splitting.  That  is,  the  image  current  appears  on  a number  of  lines 
even  though  the  image  charge  is  located  under  the  zeroth  line  only. 
Experimental  point-source  data  from  the  University  of  Rhode  Island 
is  shown  in  Figure  V-4.  These  characteristics  are  similar,  with 
respect  to  image  splitting,  but  with  regard  to  the  apparent 
assymmetry,  the  necessary  instrumentation  was  not  available  to 
know  the  actual  location  of  the  point  source  with  respect  to  the 
scan  lines.  Therefore,  these  data  may  provide  misleading  compari- 
sons regarding  assymmetry  with  the  response  of  the  model. 

Another  feature  of  assymmetry  is  shown  in  Figure  V-3.  There  is 
considerably  more  peak  current  read  out  on  the  line  before  the 
center  line  compared  to  that  on  the  line  after  the  center  line. 

With  a slower  read-out  velocity,  as  shown  in  Figure  V-5,  the  main 
peak,  in  fact,  occurred  one  line  early.  In  Figure  V-6  where  the 
read-out  speed  was  the  same  as  that  used  to  produce  Figure  V-3, 
the  point  appears  to  have  been  half  way  between  the  center  line  and 
the  preceding  line.  Actually,  the  image  charge  was  located  on  a 
diode  only  one  quarter  of  the  distance  from  the  center  line  to 
the  preceding  line.  This  feature  illustrates  a flaw  in  the  usual 
experimental  procedure  in  which  a point  source  is  centered  on  a scan 
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EXPERIMENTALLY  MEASURED  POINT  SOURCE  RESPONSE 
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Figure  V-5  POINT  RESPONSE 
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Figure  V-6  POINT  RESPONSE 


line.  In  the  procedure,  the  scan  line  locations  are  adjusted 
so  that  the  equal  height  peaks  are  obtained  on  the  two  lines 
adjacent  to  the  one  containing  the  main  peak.  The  model  shows 
that  without  background  (discussed  in  the  next  section) , this 
procedure  does  not  physically  locate  the  point  source  on  the  target 
(called  mechanical  boresighting) . 

In  all  of  the  above  cases  the  image  was  confined  to  one  diode 
(the  smallest  possible  image)  and  the  background  was  set  to  zero. 
However,  response  to  images  distributed  over  a number  of  diodes  is 
of  interest.  In  particular  the  ratio  of  the  mean  signal  current 
to  the  mean  flat-field  current  as  a function  of  image  area  on  the 
target  provides  insight  in  cases  where  images  are  blurred  by  optical, 
atmospheric  and  other  effects.  This  ratio  will  show  a maximum  image 
area  beyond  which  the  contribution  to  read-out  current  is  no  longer 
affected  by  increased  image  area.  These  data,  obtained  from  the 
Night  Vision  Laboratory  of  the  Army,  were  taken  in  which  the 
illumination  density  was  dept  constant  but  the  spot  (image)  size 
was  increased.  The  same  simulated  data  were  taken  from  the  model 
and  they  are  compared  in  Figure  V-7. 

It  is  seen  that  the  read-out  process  behaves  as  an  energy 

2 

analyzer  for  spots  of  less  than  .004  to  .008  m in  area  and  as 
an  energy  density  analyzer  for  larger  spots.  In  the  simulation, 
the  behavior  was  the  same  and  practically  the  same  spot  area  was 
obtained  for  the  transition  from  one  regime  to  the  other.  The 
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Figure  3ZI- 7 SPOT  IMAGE  RESPONSE/ BROAD  AREA  RESPONSE  vs  SPOT  SIZE 
(EXPERIMENTAL  DATA  FURNISHED  BY  R.  GRAFT,  ARMY  NIGHT  VISION  LAB) 


differences  observed  are  due  to  the  fact  that  the  read-beam 
diameter  was  not  determined  in  the  experiment;  it  was  estimated 
and  set  in  the  simulation  to  be  50  ym,  which  may  have  been  too 
small . 

5.3  POINT  SOURCES  IN  THE  PRESENCE  OF  BACKGROUND  NOISE 

The  more  realistic  situation  is  where  there  is  a point  source 
which  is  to  be  detected  but  is  observed  in  the  presence  of  back- 
ground noise.  Figure  V-8  shows  the  simulation  of  this  situation 
where  the  intensities  of  the  image  and  background  were  chosen  to 
represent  the  most  difficult  GEODSS  case.  The  point  source  con- 
sisted of  an  equivalent  of  130  electrons  released  from  the  photocathode 
to  one  diode  and  the  background  consisted  of  the  equivalent  of  an 

g 

average  of  2 x 10  electrons  released  by  the  total  photocathode  area. 
The  presence  of  background  changes  the  point  response  consider- 
ably in  a number  of  ways.  First,  it  is  seen  that  the  response  is 
more  symmetrical.  The  peak-current  location  on  the  current  map 
occurs  more  nearly  to  that  of  the  location  of  the  center  of  the 
image  charge  along  a scan  line  (Figure  V-8).  Secondly,  as  shown 
in  Figure  V-9 , when  the  point  image  is  moved  one-quarter  the  distance 
from  the  center  line  towards  the  preceding  line,  the  largest  peak 
still  occurs  on  the  center  line.  This  is  quite  a different  response 
than  that  shown  in  Figure  V-6  where  the  image  was  moved  to  the  same 
location,  -ut  where  there  was  no  background.  It  indicates  that  but 
the  point  response  is  more  symmetrical  in  the  vertical  direction  as 
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Figure  V-8  POINT  SOURCE  PLUS  BACKGROUND 


% 


r 


a 


X 

a 

i 

x 


V 

'o 


84 


Figure  V-9  POINT  SOURCE  PLUS  BACKGROUND 


well.  Figure  V-9  illustrates  this  fact  better  than  Figure  V-8, 
although  it  is  somewhat  visible  there  too,  but  is  obscured  by  the 
response  to  the  background  noise. 

Distributed  point-source  response  was  also  of  concern  for 
reasons  stated  earlier.  An  image  corresponding  to  a total  of  130 
photoelectrons  was  distributed  using  three  simulation  runs  from 
one  diode  to  four  and  then  to  nine  diodes  to  observe  the  change  in 
shape  of  peak  current  as  well  as  the  change  in  the  magnitude  of  the 
peak.  As  shown  in  Figure  V-10,  where  the  image  is  distributed 
over  four  diodes,  there  is  little  change  in  either  feature  but  it 
was  observed  that  there  was  considerable  change  when  the  image  was 
distributed  over  9 diodes. 

While  the  change  in  point  response  as  a result  of  adding 
background  is  significant,  probably  the  most  important  consider- 
ations are  those  of  signal  to  noise  ratio  (SNR)  , probability  of 

detection  (P_)  and  probability  of  false  alarm  (P.,.)  caused  by  peaks 
D f A 

in  the  background  noise.  The  ability  to  determine  these  is  crucial 
to  overall  GEODSS  system  design,  and  the  statistical  evaluation 
programs,  described  in  Section  4.2.2,  were  and  shall  be  used 
extensively  to  determine  these  numbers.  Examples  of  the  use  of 
this  program  will  illustrate  how  they  are  determined. 

As  indicated  earlier,  through  a separate  analysis  (not  published), 
the  probability  density  function  (PDF)  of  the  current  contribution 
from  a single  diode  was  computed  and  plotted.  It  was  found  that  its 
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Figure  V-10  EXTENDED  POINT  SOURCE  PLUS  BACKGROUND 


shape  matched  the  Gaussian  PDF  shape  very  closely  over  a wide  range 
of  illumination  and  read-out  speeds.  Since  the  read-out  beam 
irradiates  a number  of  diodes  causing  their  current  contribution  to 
be  weighed  and  summed,  the  PDF  of  the  sum  of  the  contributions  will 
approximate  a Gaussian  PDF  more  closely.  This  was  born  out  experi- 
mentally by  the  University  of  Rhode  Island  as  shown  in  Figure  V-ll. 

Hence,  with  the  ensemble  means  and  standard  deviations,  computed 
by  the  statistical  evaluation  programs  at  a point  where  the  peak 
image  signal  is  located,  the  equation  for  ensemble  PDF’s  can  be 
written.  Two  cases  can  be  considered;  one  in  which  only  background 
noise  is  present  and  the  other  where  an  image  is  also  present.  SNR 
can  also  be  determined  the  results  of  these  two  cases  as: 

(Mean  of  signal  + background)  - (Mean  of  background)  (5.1) 

(standard  deviation  of  signal  + background) 

The  PDF’s  of  Figure  V-12  and  V-13  were  obtained  in  the  way  just 
described  using  a set  of  typical  tube  and  operating  parameters  and 
with  two  background  illuminations  but  with  the  same  image  intensity. 

In  these  examples,  if  as  shown  in  the  figures  the  detection  thres- 
hold is  get  so  that  a P^  of  .983  required  by  the  GEODSS  system  is 

achieved,  there  is  a drastic  difference  in  the  P ’s.  In  Figure  V-12 

rA 

the  P is  .00006  whereas  in  Figure  V-13  the  P is  .0323,  a dif- 
r A r A 

ference  of  more  than  2 1/2  orders  of  magnitude.  In  both  figures, 

P^  is  equal  to  the  area  to  the  right  of  the  threshold  line  under 

the  P^  curve  on  the  right  while  P^.  is  equal  to  the  area  to  the 
Dr  r A 
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-II  MEASURED  PDF  OF  OUTPUT  NOISE 
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12  PDF'S  OF  SIGNAL  + BACKGROUND  AND  BACKGROUND  ALONE 
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Figure  3ZI-I3  PDF'S  OF  SIGNAL  + BACKGROUND  AND  BACKGROUND  ALONE 


right  of  the  threshold  line  under  the  P_„  curve  on  the  left. 

Dr 

A capability  of  the  model  is  the  ability  to  easily  relate  SNR 
of  the  output  current  back  to  the  SNR  on  the  target,  thereby 
defining  a maximum  area  on  the  target  affecting  the  output  SNR.  As 
described  in  Section  4.2.2,  the  mean  and  variance  of  the  current  is 
computed  for  each  diode  and,  since  each  contribution  is  independent, 
they  are  summed  to  get  the  total  mean  and  variance  of  the  target 
current.  Hence  each  diode  contribution  is  individually  identified 
in  the  model,  which  is  not  the  case  in  the  usual  operation  of  an 
actual  EBSICON.  Herein  lies  a powerful  capability  of  the  model 
when  exercised,  for  example,  as  described  in  the  following  investi- 
gation. 

The  authors  were  interested  to  know  whether  there  existed  a 
minimum  area  on  the  target  so  that  the  noise  from  any  larger 
area  had  no  effect  on  the  SNR  of  the  output  current.  So,  instead 
of  summing  all  the  weighted  contributions  to  determine  SNR  of  the 
output,  separate  specialized  computer  runs  were  made.  SNR  was 
calculated,  using  equation  5.1,  but  only  contributions  from 
subarrays  were  used  to  obtain  the  means  and  standard  derivations. 
First  a one  diode  array  was  used,  then  a 2 x 2 diode  subarray  was 
used  continuing  to  a 3 x 3 subarray,  etc.  Each  subarray  contains 
the  signal  peak.  In  Figures  V-14  and  V-15,  the  resulting  SNR  is 
plotted  as  a function  of  the  number  of  diodes  on  a side  of  each 
square  array.  It  is  seen  that  in  both  cases  each  curve  reaches 
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-14  SNR  TRANSFER  FROM  OUTPUT  TO  TARGET 
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-15  SNR  TRANSFER  FROM  OUTPUT  TO  TARGET 


the  final  SNR  when  the  contributions  are  summed  from  a subarray  of 

2 

5x5  diodes.  This  area  (.004  mm  ) corresponds  remarkably  well 

2 

with  that  determined  in  Figure  V-7  (.004  to  .008  mm  ).  Thus  con- 
tributions from  larger  subarrays  do  not  affect  the  SNR  of  the  out- 
put. It  is  important  to  define  the  maximum  area  on  the  target  that 
need  be  considered  for  SNR  computations.  Considering  areas  larger 
than  this  yields  a lower  SNR  on  the  target  which  would  require  a 
larger  telescope  to  regain  or  improve  the  SNR.  Further  consider  the 
implications  of  Figure  V-14 , for  example.  If  the  output  SNR  of  5.8 
is  acceptable  for  the  GEODSS  system,  then  the  required  SNR  on  the 
target  need  be  only  3.5  (not  5.8).  This  reduction  has  a large 
impact  upon  the  sizing  of  the  optical  system  used  to  illuminate  the 
image  tube.  In  this  case  the  reduction  represents  a 40%  reduction 
in  the  telescope  diameter.  The  target  area  defined  above  is  most 
likely  signal  and  background  dependent,  but  nevertheless,  the  above 
capability  can  be  used  to  easily  examine  many  anticipated  conditions 
and  which  is  not  possible  to  do  using  an  actual  tube. 

Another  feature  of  the  model  described  in  Section  4.2.2  is  the 
ability  to  obtain  the  normalized  autocovariance  function  of  the 
video  output  noise  shown  in  Figure  V-16  for  typical  tube  and  oper- 
ating parameters.  The  Fourier  transform  of  this  function  yields  the 
power  spectral  density  and  the  intrinsic  bandwidth  of  the  video 
noise . 
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Figure  2-16  NORMALIZED  AUTOCOVARIANCE  FUNCTION 


5.4  RESOLUTION 


Resolution  is  usually  measured  by  illuminating  the  photocathode 
with  a bar  chart  consisting  of  vertical  or  horizontal  bars  alter- 
nating from  solid  black  to  solid  white  and  extending  in  length  from 
near  the  bottom  to  near  the  top  of  the  picture  or  from  one  side  to 
the  other  (called  TV  lines) . A typical  output  current  map  is  shown 
as  simulated  in  Figure  V-17  for  a particular  bar  chart.  The  output 
current  is  measured  for  specific  spacings  of  the  bars  (i.e.,  the 
number  of  TV  lines  per  raster  height) . As  the  number  of  TV  lines 
per  raster  height  is  sequentially  increased  (i.e.,  different  bar 
charts  with  closer  spaced  bars),  the  maximum  minus  the  minimum  out- 
put current  is  divided  by  the  maximum  plus  the  minimum  output  current. 
The  resultant  function  is  called  the  square  wave  response  or  modulation 
transfer  function  (MTF) . 

Comparison  of  the  model  was  made  with  two  sets  of  experimental 
MTF  data  shown  in  Figure  V-18.  The  curve  marked  with  circles  shows 
Night  Vision  Laboratory  data  of  a tube  with  a similar  but  smaller 
target.  The  curves  marked  with  xfs  and  triangles  are  simulated 
MTF  data  using  a 65  ym  and  50  ym  read  beam  respectively.  The  dashed 
curve  is  the  present  estimate  of  the  new  Westinghouse  80/32  tube 
(Section  II)  MTF.  It  can  be  seen  that  the  model  has  the  appropriate 
trends  and  that  by  adjusting  parameters  in  the  model  that  it  can  be 
made  to  correspond  well  with  experimental  data. 

Resolution  would  also  include  a measure  of  the  ability  to 
distinguish  between  two  point  sources  separated  by  some  minimum 
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Figure  V-17  TYPICAL  RESPONSE  TO  A BAR  CHART 
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Figure  Z2T-I8  HORIZONTAL  SQUARE  WAVE  RESPONSE  OF  DIFFERENT  TARGETS  (EXPERIMENTAL 
DATA  FURNISHED  BY:  R.  GRAFT,  ARMY  NIGHT  VISION  LAB  AND  WEST  CORP.  HORSE  HEADS  N.Y. ) 


distance.  Currently,  the  authors  have  no  experimental  data  on  point 
source  resolution.  The  simulation  was  used  with  a set  of  typical 
tube  and  operating  parameters  to  examine  two-dimensional  resolution. 
Two  point  images  each  representing  the  equivalent  of  130  photo- 
electrons released  by  the  photocathode  were  localized  each  to  a 
diode  on  a scan  line  and  in  three  separate  simulation  runs  were 
separated  by  one,  two  and  then  three  diodes.  Figures  V-19 , V-20 
and  V-21  respectively  show  the  output  current  response  maps  for 
these  three  conditions.  As  seen,  with  the  50  ym  beam  used  in  the 
simulation,  three  diodes  of  separation  are  required  to  clearly  dis- 
tinguish the  images  along  a scan  line  in  the  presence  of  background 
noise  with  an  average  of  50  photoelectrons  per  diode  from  the  photo- 
cathode used  in  the  simulation.  Since  the  points  can  be  placed 
diagonally  or  vertically  spaced  as  well,  a more  complete  investigation 
would  be  to  include  the  Se  cases  also.  The  authors  plan  to  accom- 
plish this  when  more  actual  tube  and  operating  parameters  become 
available  for  the  new  Westinghouse  80/32  tube, 
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Figure  V-19  POINT  SOURCE  RESOLUTION  WITH  ONE  DIODE  SEPARATION 
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Figure  V-20  POINT  SOURCE  RESOLUTION  WITH  TWO  DIODE  SEPARATION 
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Figure  V-21  POINT  SOURCE  RESOLUTION  WITH  THREE  DIODE  SEPARATION 


SECTION  VI 


SUMMARY  AND  RECOMMENDATIONS 


6 . 1 SUMMARY 

The  scene  to  which  an  EBSICON  TV  camera  tube  is  exposed  depletes 
charge  stored  on  the  target  of  the  tube  in  proportion  to  the  intensity 
of  illumination  over  the  two-dimensional  scene.  As  discussed  in 
Section  III,  this  exposure  process  involves  the  statistical  phenomena 
associated  with  the  photocathode,  the  gain  mechanism,  and  the  charge 
storage  mechanism  of  the  tube  target.  After  exposure,  the  tube  is 
ready  to  be  read  out  and  the  read-out  process  provides  the  output 
video  signal  of  the  tube  as  also  described  in  Section  III.  The  read- 
out process,  as  well,  restores  the  charge  depleted  during  exposure. 

The  translation  of  the  mechanics  of  the  read-out  process  into 
closed-form  mathematical  equations  forms  the  heart  of  the  mathematical 
EBSICON  tube  model.  The  remaining  important  features  of  the  model  are 
the  mathematical  descriptions  of  the  statistical  phenomena  associated 
with  the  exposure  process. 

The  model  is  structured  on  the  basis  of  contributions  to  the 
total  output  current  signal  from  individually  identifiable  target 
diodes.  Each  contribution  is  a function  of  diode  capacitance,  tube 
gain,  average  input  signal  and  background  intensities,  read-out 
beam  width  and  current,  read-out  scan  speed,  target  raster  height  and 
width,  read-out  scan-line  density,  the  voltage  at  which  read-out  beam 
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limitation  occurs  and  the  illumination  intensity  at  which  target 
saturation  occurs.  These  are  all  adjustable  parameters  in  the  model 
thereby  allowing  the  tube  designer  to  predict,  before  building  and 
testing  a tube,  the  performance  of  the  tube  for  various  tube  and 
operating  parameter  values  and  under  various  illumination  conditions. 
Critical  performance  trade-offs  which  can  impact  cost,  can  be  made. 

In  addition,  due  to  the  structure  of  the  model,  it  can  provide 
information  not  easily  measured  (nor  usually  measureable)  in  the 

laboratory.  An  important  example  of  this  is  described  in  Section  5.3. 

In  the  model,  as  just  indicated,  each  target  diode  contribution  to  the 
total  output  current  can  be  individually  identified.  Using  this 
capability,  the  authors  were  able  to  determine  the  minimum  area  on 

the  tube  target  which  contributed  to  the  output  signal  to  noise  ratio 

(SNR) . By  knowing  this  area,  the  required  SNR  of  the  input  was  found. 
The  SNR,  in  turn,  determines  the  size  of  the  telescope  required  for  the 
system.  The  model  is  a highly  flexible  tool  useable  for  the  under- 
standing of  the  behavior  of  EBSICON's  and  for  translating  tube  and 
operating  parameters  into  Performance  information. 

As  described  in  Section  4.2,  the  model  is  implemented  in  three 
computer  programs.  One  program  simulates  the  actual  video  output 
current  signal  of  the  tube  and  as  described  in  Section  V is  useable 
in  determining  resolution  of  two  point  sources  in  the  presence  of 
background,  the  square  wave  response  of  the  tube  and  the  point  spread 
function  of  the  tube.  The  other  two  programs  are  used  to  evaluate 
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the  statistical  performance  of  EBSICONfs.  The  outputs  of  these  pro- 
grams are  used  to  determine  the  probability  of  detection  of  various 
images  in  the  presence  of  background  noise,  the  probability  of  false 
alarms  due  to  the  noise  and  signal  to  noise  ratio  for  various  back- 
ground noise  levels. 

Section  V,  in  general,  compares  the  results  obtained  from  the 
model  with  experimental  EBSICON  data  that  was  available  at  the  time. 

It  is  shown  there  that  all  trend-like  results  of  the  model  agree  with 
experimental  data.  Further,  for  each  experiment  in  which  tube  and 
operating  parameters  were  better  known,  then  when  they  were  set  to 
the  same  values  in  the  model,  quantitative  agreement  improved 
greatly.  From  this,  it  is  felt  that  the  implementations  of  the 
model  encompass  every  EBSICON  characteristic  and  operating  mode  of 
interest  anticipated  at  this  time. 

The  system  designer  can  use  the  model  by  varying  the  parameters 
in  a systematic  manner  to  achieve  predicted  optimum  overall  system 
performance.  Simultaneously,  with  the  tubes  available,  pertinent 
experiments  should  be  performed  to  insure  agreement  of  results  from 
the  model  and  from  the  experiments.  Examples  included  in  this  effort 
would  be  the  determination  of  optimal  read  beam  size  and  current, 
scan  speed  and  gain  magnitude  for  various  images  and  background  levels 
to  achieve  maximum  probability  of  detection  and  minimum  probability 
of  false  alarms  while  at  the  same  time  retaining  compatibility  with 
reasonable  data  rates.  Having  arrived  at  an  optimal  set  of  parameters, 
the  required  tube  and  camera  system  can  be  specified. 
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6 . 2 RE  COKMENDAT  TONS 


Currently,  under  the  GEODSS  program,  the  Night  Vision  Laboratory 
(NVL)  of  the  Army  is  conducting  a tube  development  program.  Vest- 
inghouse  in  coniunction  with  Texas  Instruments  is  providing  EBSICON 
image  tubes  and  yokes  to  NVL  to  be  tested  using  NVL’s  laboratory 
camera  system.  The  experiments  will  provide  measured  tube  character- 
istics. Manv  characteristics  or  parameters  such  as  beam  width, 
target  saturation  illumination  and  the  illumination  value  causing 
read-out  beam  limitation  cannot  be  directly  measured  with  an  actual 
EBSTCON.  In  fact  the  effects  of  the  last  two  characteristics  are 
difficult  to  separate  by  measurement.  Many  characteristics  must 
be  estimated  which  leads,  in  many  cases,  to  inaccurate  values  of 
the  actual  physical  parameters  of  the  tubes.  The  tube  model  can 
be  used  to  support  in  improving  the  accuracy  of  the  tube  parameter 

estimates . 

Operating  parameters  such  as  read-out  speed,  tube  gain  and  read- 
out beam  current  determine  how  the  tube  and  camera  system  effect  the 
overall  system  performance.  The  influence  of  these  parameters  on 
system  performance  can  be  estimated  from  the  model  and  a set  of 
recommended  operating  settings  can  be  determined. 

Finally,  the  model  can  be  used  to  assist  in  formulating  the 
acceptance  criteria  for  the  image  tube  and  camera  systems  which  will 
be  used.  Specifically,  it  is  recommended  that  the  image  tube  model 
be  used  to: 
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1.  Support  the  NVL  activity.  This  support  will  consist  of 
two  activities. 

First  the  model  will  be  used  to  verify  those  tube  parameters  which 
can  only  be  estimated.  This  shall  be  done  by  setting  the  parameters 
in  the  model  to  the  estimated  values  for  a set  of  experiments  which 
have  been  or  will  be  performed  that  clearly  demonstrate  the  effects 
of  the  actual  values.  The  operating  parameters  in  the  model  are  set 
to  the  values  used  in  the  experiments,  and. the  model  is  then  exercised. 
Since  these  parameters  are  clearly  identified  in  the  model,  they 
will  be  varied  systematically  in  the  model  until  agreement  of  results 
from  the  model  and  from  the  experiments  is  good.  The  resultant  values 
obtained  from  the  model  will  provide  accurate  estimates  of  the  actual 
values  of  the  physical  parameters  of  the  tube. 

Second,  using  the  model  as  an  analytical  tool,  additional 
experiments  designed  to  improve  the  evaluation  of  the  tube  parameters 
will  be  identified  and  recommended. 

2.  Specify  Operating  Parameters.  For  tubes  with  acceptable 
physical  parameters,  the  model  will  be  used  to  establish  the  optimum 
set  of  operating  parameters  (such  as  read  speed,  gain,  etc.)  which 
yield  the  required  system  performance.  In  all  cases  performance 
requirements  will  be  based  upon  probability  of  detection,  probability 
of  false  alarms,  lag,  resolution,  blooming,  etc.  and  shall  be  com- 
patable  with  reasonable  video  signal  bandwidths.  The  model  will  be 
used  to  establish  these  factors,  but  they  shall  also  be  shown  experi- 
mentally for  each  tube  under  selected  illumination  conditions. 
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3.  Generate  Acceptance  Criteria  and  Performance  Requirements 
for  Image  Tube  and  Camera  Systems.  The  interrelation  between  the 
tube  model  and  the  tube  development  program,  as  discussed  above, 
shall  be  the  basis  on  which  the  EBSICON  tube  and  camera  system 
acceptance  criteria  are  generated  as  part  of  the  overall  system 
specification.  From  the  knowledge  of  EBSICON  performance  and  its 
impact  upon  the  overall  system,  performance  requirements  for  other 
types  of  tube  and  camera  systems  will  be  generated.  These  will  be 
of  a generic  nature  such  as  point  source  resolution,  modulation  trans- 
fer function,  two-dimensional  point  spread  function,  lag,  dark  current, 
sensitivity,  microphonic  effects,  and  signal  and  noise  transfer 
relationships.  All  of  these  can  be  determined  by  laboratory  measure- 
ments which  are  strongly  influenced  by  experience  with  the  tube  model 
and  EBSICON  measurements. 
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APPENDIX  A 


SIMULATION  PROGRAM  LISTING 
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./  ADD  SEQFLD=738 

./  NUMBER  NEW  1 -00000001 , I NCR =00 0000 01 

REAL»8  U/. 3 00 3 00 3003/ _ 

IEB801I  NM  BLOCKSIZE  ASSUMED  80. 

DIMENSION  A(11,101),Q(21,21),Q0(21,21) 

INTEGER  LN(3)/6,8, 10/ 

COMMON  IBUF(IOOO) 

REAL  NL(3 )/600. ,800.  , 1 000  . / , L , 10 , LN2 

REAL  XT ( 3 ) / * X-AX  * , »IS  ft  »CH>  • / , YT  ( 2 W • L INEJ  , • NO.  * / , Z T ( 4 ) / • CURR  1 

1 1 1 ENT  1 , • ( AMP  • t • S ) •/ 

__  REAL  JNE(2)/1.6019E-14f  8.0095E-I5/  __  _ 

VMB  = 1 9 • 

VMT  =16. 

G A IN= 1 00  . 

NET  = . 32*10.  ♦ *(9.-.<»*VMT) 

10=1. 0E-6 

ALPHA=.0025  _ 

BETA=‘.25£+15 

P 1 = 3 .141  59265 

PIR  = SORT (PI) 

DA  = 9 . 37  5E“7 

FAC J=IO*DA / ( P I * ALPHA  **2 ) 

H =2 . 5 
T=l/17.5 

W =H  

X0=-.025 
1 = 1 

V=NL ( t ) *W/T 
TS  = #02  *T /NL ( I ) 

L = H / N L ( I ) 

LNH  =LN ( I ) /2  . 

LN2  =LNH 
E PS  =0 .0 
DO  70  J = l , 2 

DO  5 11=1,21 __  _____ 

DO  5 JQ  = 1 ,21 

Q( 10, JU) =0,0 

5 CONTINUE 
DO  70  K = 1 , h 
DO  10  IQ=l ,21 
1016  DO  10  JQ  = 1 ,21 

Q ( I Q , JQ ) =Q ( I Q , J qT+N E ( J ) 

I F ( K . NE . 6 ) GO  TO  10 

0(10, JO)=NE(J) 

10  QQ ( IQ , JO ) =Q ( IQ, JQ) 

WRITE (6, 1011) 

1011  FORMA  T ( __  • _ __ 

DO  1010  J0=  1 , 21 

WRITE (6, 1071)  (Q( IQ, JQ),IQ=1 ,1 M __ 

WRITE ( 6, 1072 ) (0( 10, JO ,10=1? ,21 ) 

JlOOp  CONTINUE 

LNI=LN( I )♦! 
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LNI J =2 1 *LN ( H + 1 

DO  60  N-l.LNU 

Y0=  ( M-LN ( I )-l  )*L»£PS _ 

Nl-N-LNH 

DO  40  IT=1 ,201 

TIME=(IT-1 ) *T  S/200 . 

IM=IT-50 

X=V*TTMC+X0 

C=0.0 

IY0=1 1 .5+/0/. 00125 

1X0  = 11. 5+X/. 00125 

DO  35  10-1,21 

ISQ-IQ-IX3 

I F ( IABS( ISO) .GT.S )G0  TO  36 
XD«CIO-n  )*. 0012 5 

XI-XO-X 

XIO^XD-XO 

F=-( XI/ ALPHA) **2 

FAC  X = 0 .0 

IFtE.GT. (-30. ) )FACX  = cXP(E  ) 

DO  JO  JO-1  ,21 

IF~((  ( J3-IY0)**2+ISU**I  ) .or. 66  ) iC  TO  30 
YD=( JO-1 1 ) *.001 25 
YJ= YO-YO 

E=— ( YJ/AL°HA ) **2 

FAC  Y=0.0 

IF(E.5T.  (-pn.  ) ) FAC Y=  EXPCI  ) 

I F ( 0 ( 10, JJ)  ) 3 0 * 3 0 , 2 0 

20_  D0=  ( I 0*cACY*r<  A/(2.*V*ALP.1A*PIR  ) ) * ( ERF ( XI  /ALPHA  ) — trtF  ( A IO/ALPHA  ) ) 
FAC  1 — 1 .o 

Q ( I J , J >J  ) = A L 00  ( E X P ( R F.  TA  * ( 3 0 ( I J , J 0 ) 03)  )-lX  P(QTA»DJ)  + 1 . I/BCTA 
C=C+(1  .-EXP (-QSTA *3 ( I 0,J3)  ) )*rACJ*FACI*FACX*FACY 
30  CONTINUE 
35  CONTINUE 

I F ( IM.LT.l .OK. IM.GT. 101 .JK.N1 .LT.1.IK.N1 .GT.LNI ) 00  TO  60 

A (Ml, IM)=C 

60  CONTINUE 1 

DO  50  I3-],21 

DO  60  JO-1,21 

50  Q0( IQ, J0)=0( 10, JO) 

60  CONTINUE 

WRI  Tc (6,61  ) J 

61  FQRMATt * 1 J = » , 12/) 

DO  1060  IT-1,101 

1060  WRITE  (6,1061  HAINI  , IT)  ,N1  =1  ,LN() 

1061  FORMA  T ( * • ,1P51 1 ,6,10( IX, 1PE1 1.6)) 

WRITC (6, 1062) 

1062  FORMATt 1H1  ) 

DO  1070  J '3  = 1 , 2 1 

WRITE (6, 1071 ) (0(1 j, JO) ,1 3 = 1 ,1 1 ) 

1070  WRITE (6, 1072) (0(10, JO), I 0-1 2, 21) 

1071  FORMATt*  » , 1 1 ( l PE  1 o . 3 ) ) 
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1072  FORMAT! 1 ' . IPX , 10 (I P E10 . 3 ) ) 

70  CONTINUE 

CALL  EXIT 

END 

FUNCTION  RAND ( X ) 

REAL*?  Yf X ' 

DATA  X*V 10000  3 ,/ 

Y=X*XM 

M=Y _ 

X=Y-M 

RAND=X 

RETURN 

END 

FUNCTION  GAUS(Ul,U2iVMB) 

AM»2.*10.**19.-.4*VMB) 

" GAUS=SQRT  ( -2  . *AM*\LOG(  U.l ) ) *COS<  6.2B195*J2  H-AM 
RETURN 
END 

IEB818I  HIGHEST  CONDITION  CODE  WAS  00000000 

IEB819I  END  OF  JOB  IErtUPDTE. 
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APPENDIX  B 


STATISTICAL  EVALUATION  PROGRAM  (NO. 


LISTING 
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./  ADD  SEQF  LD  = 738 

./  NUMBER  NEW1 =00000001 , I NCR =00000001 

REAL*9  AltA2t  A3fBl»B2»B3»EIl»EI  12, El  1_1 2 , COV 

IEB801I  NM  BL0CKSI2E  ASSUMED  80. 

1,PIR2,FACX1 ,FACX2 ,FACX3,DQ,  GI R , E PS  , FMAX  ,_VMB  _ 

REAL*8  ARG,ARG1 

REAL*8  FACY,IBlti)Ql»Klt M3 , M4,M5 , M6 , F AC X3l , F ACX21 

REAL* 8 IOtNLt LtIB,Kt II , 12 , 1 R , MU , ME AN , AM, PI,FAC1 , F AC 2 t S I G , V , VAR , X , Y 
REAL*8  VMT  , QT  ♦ MUT » S I GT » Ml T,M2T,M3T ,M4T , M5T, M6T tNET 
REAL*8  EXtEXE tHKGfHNGtHKGl ,H.NG1 , F,F2,F1 , FI  2 

REAL*e  ZltYNyYNltZ 

REAL*8  A,GAM'A,AR 

RE A L*£  M1(20,13),G(20,1?) ,M2(20,13) 

REA  L*8  DOQ,KQ, F0,FQ2 , AKG2 , ARG3 , HOG ♦ HUG2 , M3 ( 20, 13 ) 

DALOG  ( Z ) =2 ♦ 77* ( I X TUV  ( Z , 1, 7) -64.  ) 

GAIN=500. 

VMB  =1 9 * 

A1=2.*10.**(9.-.4*VMB) 

81=GAIN*1 .60190-19 

ALPHA=.0025 
BET  A =3 • E 1 3 
G A M A = 4 • 4 2 8 5 8 E + 1 1 
A = 8 ET  A/G AM  A 
X0=-.025 

P 1=3 . 141 532  65  3 59 9 7 93  POO 

D A=  9 . 37  S E -7 
1 0 = ? * F — 6 

FAC1=I0*DA/ ( PI* ALPHA **2 ) 

H=?_.5 

NL=~1000. 

W=H _ 

N0=0 

L=H/NL 

PIR?=DSQRT(?.*PI ) 

T=1 ./30. 

V=NL*tf/T 

FAC?=I0»DA/(2.»V»ALPHA»DSQKT(  PI  ) ) __ _ __  _ _ 

I ST0P=1 3+NO 
DO  3000  I I STOP 
DO  3000  J = 1,13 
M2(  UJ)  =0.0 
M3 ( I , J)=0.0  ' 

3000  CONTI NUC 

00  3010  1 = 1 , 1 STOP 
DO  3010_J=1,13 
Ml ( I t J ) =A1 
G ( I * J ) =8 1 

MU=( 1 .-DEXP (-Ml ( I , J ) *(  1 .—DEXP (-GAMA*B1 ) ) ) )/GAMA 

VAR=( DEXP (-MI  ( I , J ) * ( 1 .-DEXP  (-2.  *GAMA*B1  ) ) ) -D.3  XP  ( -2  . *M  1 ( I , J ) * ( 1 . -Dr 

1 XP ( -G AM A*B !))))/( GAM A** 2 ) 

VAR  =V  AR  + M3 ( I,  J ) 

MU=MU+M2  ( I , J ) 
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Ml  ( I , J ) = (0L0G(1.-GAMA»MU)**2)/ULUG(( l.-GAMA«MU)/(VAR*GAMA«*2♦(l.- 
l GAMA*MU ) **2 ) ) 

G ( I , J ) =-PLOG ( DLOG ( VAR*GAMA**2  + ( 1 ,-GAMA*MU ) **2 )/DLUG(  1 ,-GAMA*MU)-1  ♦ 
D/GAMA 

3010  CONTINUE  _ 

MEAN=0.0 

VAR =0,0 

C0V=0.0 

WRITE (6, It ) 

11  FORMAT! • 1 • » 8X» 'K*  1 13X* ’03* i 9X,'M122',  9X , • M 1 2 • , 1 OX , * M2 • , 10X , 'Ml  • , 

1 tOX,' ME AN', IPX, 'VAR* ,1QX,  »C0V',7X,'?  TERMS  LOST') 

Ob  100  I =1 « I STOP 

X = ( 1-7)*, 00)25 

FACX1 =FAC7*(DERF( ( X+XO ) /A LPHA l-DERP ( ( X-XO ) / A L PH A ) ) 

FACX2=FAC2*(0ERF( X/ALPHA) — DERF ( (X-XO) /ALPHA ) ) 

FACX21=FACZ*(DFRF ( ( X-NO* . 00 1 2 5 ) / ALPHA ) -DERF ( ( X-X0-N0*,001  25 ) /ALPHA 

1)1 

FACX3=FACl*DEXP(-( X/ALPHA )**? ) 

FACX31 =FAC1*DEXP(— ( ( X-NO* ,00125 ) /ALPHA ) **2 ) 

DO  90  J=1 ,13 

Y=( J-7)*.nr>125 

FAC Y=DEXP( -(Y/ALPHA ) **2 ) 

I5t=FACX31*FACY 
IB=FACX3*F  ACY 

003=0  ,0 

DQ=FACX2*FACY 
DOT  =F AS  X2 1 «FaCY 
DO  10  N=l, 5 

DQ1=DJ1  +FACXl*DEXP(-(  ( Y+\I*L  ) /ALPHA  ) ) 

D0=DQ+FACX1*DEXP(— ( ( Y+N*L ) /ALPHA ) **2 ) 

10  CONTINUE 

DO  20  N= 1,10 

D3Q=D31*FACXl*DEXP(-( ( Y-( N— 5 ) *L ) /ALPHA ) **2 ) 

20  CONTINUE 

KO=DFXP(BETA*DQO ) 

K=OEXP(35TA*DQ) 

K 1 =DExP ( BETA*D31  ) 

AM=M1  ( I , J ) 

J.  M^AM+A.*D.rQRT(AM)  + .5 
EXE=0,0 

IFJAM.GT. ( 173 . ) ) E XE = 1 73 .-AM 
IF (M.GT.346) WRITE (5,6) 

X FORMAT ('  MEAN  PHOTOELECTRONS  GREATER  THAN  ALLOWED ' ) 

NSTOP=M 

IF(M.OT.173)NST0P=173 
EX=— A M 

IF(  FxE.LT . (o.o) )EX=-173. 

Z1=DCXP ( EXE  ) 

£ » 0 

F02=0 .0 
c=0.0 
F2=o.n 


115 


F 1=0*0 

F12=0.0 

KCHJNT  =0 

YN=DtXP ( EX  ) 

DO  2035  N=  1 »MSTOP  _ _ 

YN=YN*( AM/N) 

AR=-N»G(  It J)*GAMA 

ARG=-A 

IF ( AR.GT»(-173. ) ) ARG=-A*(  1 . -DE  XP ( AR  ) ) 

ARG1 =0.0 

IF  ( ARG.OT.  (-173.  ) ) AR  G1  =l)E  XP  (ARG) 

ARG2=1.0 

I F ( (DALQG (Kg) ♦PALOG( AKG1-1J ).GT>I-173.)  ) ARG2=KQ* ( ARG  1-1 . ) + l . 
ARG3=DL0G( ARG2 ) 

HQG=0  >0 

I F ( ( DALOG( ARG3)+DALQG( 1 . / BE TA ) ) . GT . ( -1 73. ) ) H0G=ARG3 / BETA 
HQG2=0.0 

I F ( ( 2 •♦DALOGI HOG) ) .GT. ( -173. ) )HQG2=H0G**2 

HKG=(K*( 1 . -ARG 1 ) ) / ( K » (1  .-K  ) *A RG 3 ) 

HKG1=(K1*(  1 .-ARG  1 ))/(Kl+(1.-KU  * A R G 1 ) 

HNG  =0.0 

I F ( ( 2 .*DALOG  ( HivG  ) ) .GT.  (-173.)  )HNG=iKG**Z 
HNG 1=0.0 

“ I F ( (DALOG ( HKG ) ♦DALOG (UKG1  ) ) . GT. ( - 1 7 3 . ) )HNG1=H*G*H<G1 

2003  FM=0 . 0 

XTE  $T=DALDG ( YN ) +DALOG( 1 1 ) 

I F ( XTjrST.  VT.  (-173  . ) ) FM=YN*Z1 

I F ( ( XTE  S T + DA  LUG  ( HOG  ) ).GT.  (-17  3.  ) ) FwJ-FW  + FM^H  JG 
IF  ( ( X TES  T +_D  A LOG  (H0G2  ) ) • G T . ( - 1 7 3 . ) ) FQ2  =r02  +F^*HQG2 
I«=(  ( X T E S f ♦ DALOG  (HKG  ) ).GT.  (-173.  ) )F  = F + FM*HKG 

I^(  ( XTES  r + Q ALQG  ( HKG1  ) ) .GT  . ( -1  73  . ) ) F 1 =F 1 1 

IF(HNG.jJ. (0.0) )GG  TO  1010 

I F ( ( XTrST  + OAlOG  (HNG  ) ).GT.  (-173.  ) ) F 2 =F 2 +FM *HNG 

1010  IF(HNG1 .EO. (0.0) )C0  Tn  7035 

I F ( ( XTcST+PALOG(HNGl  ) ) • GM  -1 73  . ) ) F 1 2 =F  1 2 +FM*HNGJ 

EXPON=1.20412*IABS(  I XTUVI  f iVl  i 7 ) -I  XTJV  ( FM  1 1 1 7 ) -I  XTU V ( Hi>iG f 1 t 7 j +64  ) 

IF(PXPON.LT.(  16.)  )GQ  TO  2035 

K0UNT=K0UNT+1 
* WRI Tc ( 6 »2009 ) EX PON 

200?  FORMAT ( f EXPONENT  6 I FFERt^CE**  1 tlPE13.6) 

2035  CONTINUE  _ 

I F ( M . L F • 173 )GU  TO  2090 

YN1 =1  . 

DO  20  30  ?\I=174,M 
YN1 = Y N 1 *( A M / N ) 

AR=-N*G ( Ir J ) ♦GAMA 
ARG=-A 

IF(  AR.GT. (-173. ) ) ARG=-A*(  l.-DEXP ( AR ) ) 

ARG  1 =0 . 0 

I F ( ARG.OT.  (-173.))  AR  G1  =D'rXP(  ARG) 

ARG2  = 1 .0 

I F ( (DALOGC  KQ)  ♦DALOG  ( ARG  1-1  . ) ) . GT  . ( -1  73  . ) )ARG2=»\0*(ARG1-1  . ) + l . 
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ARG3=DLOGtARG2  1 

HQG=0 • 0 

I F ( (DALOG( ARG3)+DALUG(1./BETA ) ) .GT. (-173. ) ) HUG=ARG3 /BETA 

HQG2-0. 0 

IF( ( 2 . *DALOG( HQG ) ) .GT. (-173.) )HQG2=H0G**2 

HKGS ( K* ( 1 • -ARG 1 ) ) / (K ♦( 1 .-K ) *ARG1 ) 

HKG 1 = ( K 1 *( 1 .-ARG1  ) ) / (Kl + ( 1 .-Kl ) »ARG1 ) 

HNG=0.0 

IF  C (2.*PAL0G(HKG)  ) .GT. (-173. ) ) rtNG=HKG**2 

HNG 1=0.0 

_ IF ( ( D ALQG ( HKG ) ♦DA LOG  (HKG1  ) ) .GT.  (-173.  ) )HNG1  =HKG*HKG I 
2040  FM=0.0 

XT£ST=XT£ST+DALOG( YN1  ) 

I F ( ( XT£ST+DALOG( HQG ) ).GT. (-173. ) ) FO=FO* FM*HQG 

I F ( ( XTS  ST  + PALOG (H^G2  ) ) . GT . ( -1 73 . ) ) FQ2 =FQ2+FM*hQG2 

IF(XTcST.GT. (-173. ) ) FM=YN*YN1*Z1 

I F ( (XT2ST  + DAL0G(rlKG)  l.GT.  (-173.  ) )F«Mrl^*HKG 

I F ( (XTCST  + DAL0C(HKG1  ) ) .GT . { -I  73 . ) ) FI - FI +FM*HKG1 

IF(HNG.FQ. (0.0) )GU  TO  1020 

IF ( (XTEST*DALOG(HNG) ) .GT. (-173. ) ) F2 = f=2 +FM*HNG 
1020  IF(HNG1.EQ.(Q.0))G0  TO  c03 0 

IF( (XTEST>DAL0G(HNG1  n .GT. (-1 73. ) ) FI 2-F1 2 +FiM*HNG 1 

EXP0N  = 1 . 20412*1  Aft  5 ( I XT’JV(  F2 , 1 , 7 )-I  XTUV  ( F.M  f 1 , 7 )-I  XTUV  ( HNG,  1 t7)+S4) 
IF( FXPON.LT. ( 16. ) )GU  TO  2080 

K QU  N T =K  0 U N T + I 

WRITE (6,2009) EXPON 
2080  CONTINUE 
2090  £11=0.0 

PERCNT=K0UNT*lO0. /M 

Ic( (DALOG( IB) +DAL00 ( F ) ) .GT.C-173.) )cII*IB*F 

£112=0.0 

I F ( ( 2 . *DALOC  ( I B ) ♦ DA  L 00  ( F2  ) ) • 6 T . (-17  3.  ) )EI12=Ib*lK*F2 
E 1 2 = 0 .0  _ 

IF  ( ( D A L 0 G ( IB1  )OALJG(Fl  ) ) .GT.  (-173.  ) )'  t:I2=IRI*FI 
-HI2=o_.0_  _ 

IF  ( ( DALOG  ( IB)  ♦DA LUG  ( IK  I ) +OALOO ( F 1 2 ) ) .GT.-  (-1  73.  ) ) EH  I2=I8*IB1  *FI* 

M2( It J)=F^ 

M3( 1 1 J ) =FQ2 

J. I F ( ( 2 . -DALOG  ( FQ  ) ) .GT  . ( -17  3 . ) ) M3  ( I » J ) =iM3  ( I,  J )-F3**2 

ME AN=ME AN+  2 II 

VAR  =V AR ♦ £ I T 2 

C0V=C0ViEIlI2 

I F ( (2.»DAL0G(CI1  ) 1 .G T. ( -i 73. ) ) VAR =VAR-E II **2 

IF ( ( D ALQG ( £11  )+nAL0G(EI2) ).GT.(-I73.) ) CUV -CUV -El  I *2 12 
WRITE (6, 12 )K,DQ,FI II2,u I1?»EI2,FI1 , MF AN f VAR f CUV, P> RCNT 
12  FORMAT ( • 1 ,9( 1PD13.6  ),3X,F6.3 ) 

90  CONTINUE 
TOO  CONTINUE 

STDV  = PSORT (VAR  ) 

C OR  E L =1 .0 

IF ( VAR . GT • (10.** (-70)) ) COKE L=COV /VAR 
WRITS  (6V  101  )AM,MFANI,  STDV, COREL 
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101  FORMAT  ( 1 AM  = 1 , 1 PD  1 3 . * » ' '1 E AN=  ' « 1 PD1  3 . to  » 1 ST.  DE  V . ='  . 1 PE  1 3 .6  . ' CORR 

1FLATIUN  C0EF.=',lPfl3.6) 

TIP  CONTINUE  _ 

CALL  EXIT 
END 

I E B 8 1 0 1 HIGHEST  CONDITION  CODE  WAS  00000000 

IEB819I  END  OF  JOB  I EriUPDTE  . 
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APPENDIX  C 


STATISTICAL  EVALUATION  PROGRAM  (NO. 


LISTING 
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. / ADO  ScQPLO  = 73e 

./  NUMBER  NEWt  =00000001  , INCR=OODOOr>Ol 

REAl*3  41,4^1  A3, 01  *00*33*  . II, 112, Ell  1 2,  CUV 
IEB801I  NM  BLOCK  f> I ZF  ASEUMCO  CO. 

REAL**  AR  . , AX  U 1 

REALMS  FACY,ID1 ,DU1 ,K1 , HO  , .14,  .1?  , 16  , rACX.Ol  , F AC  X2 1 

REALMS  I0,NL.L,IB,n,I1,  1]_j  IK, MJ, MEAN,  AH,  PI,  FAC  1 ,FaC2,EI  ,,V,VAR,X,  Y 

1,PIR2,FACX1 ,FACX3,FACX3,D2,  G I R » c P , F M A X , VM  B 
RE A L*f  VMT  , QT  , HUT  » SliT,  Ml  T , M.;T,M  IT  , M4T  , MBT  , "ST  , NE  T 
R = A L * P.  EX  , EXE  , HKO  , UNO,  MR-',  1 , UNCI  , F*F2,F1,C12 
REAL**  Z 1 , r N , r h 1 , Z 
REAL*)?  \,GAMA,AR 

REAL*.-  * 1 ( I ? , 1 3 ) , G M A , 1 Z ) , M2  (1  3 , 1 3 ) _ 

rc al*  ^ d.j  , ,<v , f:j, f j_  , ar,;e  , ar,,a,4  jo, m jo.:  , *3  ( 1 .» , i 5 1 
DALUG ( Z ) .77*1  TXTUV  (Z  , 1 , 7)  —.04.  ) 

GA  IN=BOO. 

VMT  = 1 6 • 

'/MB  =1 0 . 

Al  =2  . »1  0 .*»  (_o -.  4* VM  ' , | 

A 2=  . 02*10.  **  ( <i.  - . a*  V a r ) 

»1=GA IN* 1.601 90- 1 9 
DO  3000  1=1,1? 

.00  3000  J = 1 , 1 3 
1? (I , J ) =0.0 

M3( I , J) =o.O 

3000  CONTINUE-- 

AL°HA=  .OO.'S 
BETA  = '»."1'> 

GAM  A=a.  \~c'  > 3 1 - ♦ 1 1 

a=b:eta/gv-ia 

xo=- . 02. 5_ 

PI  = ?.  . 1 A 1 53  ..ho  Vs  “’  J'?1-) VOiin 

DA  = 9 . 27  B r -7 
10=2 .2-4 

FAC’=10*D\/(PI*ALP  !A**.  ) 

NL=1 000 . 

H=  ? . 5 

W=H 

ND - L = 7 
L=H/NL 

PTR2  = DC0r<T  ( 2 . * P J ) 

T=l  . / 3 0 . 

on  )oi o i=i,i a 

DO  00  10  J=i  ,13 
Ml  C I , J 1 = ■'  ’ 

I F ( I . ED.  7.  AND.. I.1-"  . 7 ) HI  ( I . J ) = M + A - 
G( I , J 1=00 

HU=(  1 .-DEXPC-M1  ( I , J ) *(  1 .-U  XP(- .,AiMA*sl  )|)|/ , A M A 

V A R = ( OeXIM-Ml  ( I , J 1 * ( \P  ( — 2 • *T  A 1 X * > 1 ) ) ) — D API—  . * H I ( I , J ) * ( 1 . - 

1XP(-Gama*B1 ) I ) j/(GAM«**2j 
MU=  M J +•  M7  ( I , J) 

V AR  = / AR  + v 0 ( I , J) 
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Ml  ( I t J ) f ( DLOo(  1 .-GAMA»MU )/DLUG(  ( 1 . -G AM A*MJ ) / ( V AR*GA MA**C  + < 1 

1GAMA*MU)**2) ) 

G ( I , J L=rDLOGiDLDG ( VAR^GAM + (1 .-GAMA*MU )**C ) /DLOG ( 1 .-GAMA*MU) -1 . 
1 ) /GAMA 
3010  CONTINUE 
ME AN =0.0 

VAR  =0 . 0 

V=NL*W/T 

FAC2=IO*DA/(2  • *V*A  LPriA*D$  J)RT ( PI  ) ) 

WRI TE ( 6, 11  ) 

1 1 FORMAT  < 1 1 1 ,PXf  • K»,  1 »0Q  1 f I OX,  fM12  1 , 12  X,  fMl ' , 12X  , •MEAN1  , 12  X , 'VAR  ' 

1 ,6X,  • S TER  M$  EXCLUDED' ) 

DO  lQQ  1 = 1,13 

X = ( I-7)*.00125 

FAC  XI  =FACr«(DERr  ( (X+XO)  /\  LPnA  ) — DERF  ( ( X-XO  ) /ALPHA  ) ) 

F ACX7=F AC ?*<DFRr(X/ ALPHA) -07 RF( (X-XO) /ALP^A ) ) 

FAC  X3  = F AC 1 *DE  XP ( - ( X/ ALPHA ) **: ) 

DO  90  J = 1 ,13 

Y=( J-7)». 00125 

FACY=DEXP  (-(  Y/ALPH  \ ) **.; ) 

Itf=FACX7*r ACY 
000=0.0 
0Q=CACX2  *F  ACY 
DO  l'n  N = 1 , 5 

oo=du+facxi*lv:xp(-(  iy+n*i ) / a l ph a ) **2 ) 

10  CONTINUE 

DO  20  N=! , 10 

Di)0=n^0fFACxi«Drxp  ( - ( ( y-(  \-o  j *l  i / alpha i**?i 
20  CONTINUE 

K0=r)“ XP  ( M!TA*DwU) 

K =D  E X p < 3 r T A »[) U ) _ 

A.M=M1  ( I t J ) 

M = A M + A. *OSURT (AM )♦. 5 
E XE  =0  . ^ 

IP(  AM.Gr.  (173.  n?X-*17C.-AM 
Ic(M.GT.?A5NRTTE(6,  M 

6 FOR  M A T ( 1 MEAN  PH  OTOE  L EC  T<  J\  c GREAT'.  < ID  V\  ALLOx-J1) 

nstop=m  

. I p ( M . G T . 1 7 3 ) N S T 0 P = 1 7 ' > 

EX=-AM 

IF(CXC.LT. (0.0) )EX=-17E. 

Z 1=D^XP(CXE ) 

FU=0 . 0 

FQ?  =0  * 0 “ 

r=0 .0 

F2  = 0 • 0 
FI =0.0 
F12=n.n 

XUUNT  = 0 _ _ 

YN”Or  X P ( u X ) 

DO  20*5  f\iT'lP 

YN=YN* ( AM/N ) 
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AR=-N*G( It J)*GAMA 

ARG  =-A 

IF( AR.GT. (-173. ) ) ARG=-A*( 1 .-DEXP( AR)  ) 

ARG1=0.0 

IF ( ARG  «GT .(—173.) ) ARGt =DS XP ( ARG ) 

ARG2=1 .0 

I F ( ( DALOG (KQ)  +DALOG  ( ARG  1-1.  ) ) .GT  . ( -173  . ) )ARG2=K0*  (AKGI-1.  )■*■!. 

ARG3=DL0G ( ARG2 ) 

HQG=0 .0  _ _ _ 

IF ( (DAL0G(ARG3)+0AL0G(1 . /BETA ) ) » GT.( -173. ) )HQG=ARG3/BETA 

HQG2=n.O 

I F ( ( 2 . *DALOG(  HQG ) ) .G~T^  (-173.  ) )H0G2«HliG**2 

HKG= ( K* ( 1,-ARGl ) ) / ( K ♦ ( 1 ,-K ) *ARG 1 > 

HNG=0.0 

I F ( ( 2 »*DALQG ( HKG ) » .GT. (-1 73. ) )HNG=HKG**2  

2003  FM=0.0 

XTEST=OALQG(YN) +DALOG( 21 ) 

IF ( ( X TEST ♦DA LOG (HOG ) ) .GT. (-173. ) ) FQ=FU+FM*HQG 

IF  ( (XTEST+DAU3G(HOG2  ) ) .GT.  (-173.  ) ) FQ2 =FQ2 +FM*hQG2 

I F ( XTEST  .GT  . (-173.  ) ) FR=YN*Z1 

I F ( (XTEST*OALQG(HKG)  ) .GT. (—17^. ) )F  = F«-FM*HKG 

IF(HNG.EQ.  (0.0)  I GO  TO  203 5 

I F ( (XTEST  + DALOGMNG)  J.GT.  ( — 1 7 A . ) ) F2  = F2 +PM*H,V, 

PXPGN  = 1 .20412*1  A6S  ( I XT'JV(  F2, 1 , 7 ) -I  XT'JV  ( FM  , T , 7 ) - I XTlJV  ( , 1 , 7 ) +o4) 

IF( EXPOV.LT. ( 16.)  )G0  TO  2033 

KOUNT  =KOUNT  + 1 
WRITE ( 6 1 2009 ) EX  PON 

2009  FORMA T ( ' EXPONENT  DIFFERENCE*' * 1 P F 1 3 . o ) 

2035  CONTINUE 

I F ( M . L p • 1 7 3 ) G 0 TO  2040 

YN1  =1  . 

DO  2030  N= 1 74  * M 

YN1=YN1*(AM/N)  _ 

AR=— N*G ( I, J)*GAMA 

ARG=— A 

IF ( AR.GT.  (-173.  ) ) ARG=-A*< "i  .-DEXP(AR  l ) 

ARG1 =0.0 

IF(ARG.6T. (-173.) JARGI =-0"  XP( ARG I 
ARG2=1 .0 

I F ( (DALOG  ( KO  ) ♦DALOG  ( ARG  1-1 .“)  ) . GT . ( -173  . ) ) AR  G2  = < 2*  (ARG  1-1  . ) ♦ 1 . 
ARG3=DL0G ( ARG? ) 

HQG=0.0 

Ic(  (OALOGt  ARG  3 ) ♦DAL0G( 1 ./&LTA)).GT.(-173.)  ) H .2G*A  RG3  / He  T A 

HQG2=0.0 

I F ( ( 2 . *DALOG  ( HUG  ) ) .GT.  (-173.  ) ) h:jG2=HQG**2 
HKG= ( K* ( 1.-ARG1) ) / ( K ♦ ( 1 .— K ) *ARG1 ) 

HNG=0 .0 

I F ( ( 2~.  * 0 A LOG  ( HKG ) ) .GT.  (-173.)  )HNG=HKG**2 

2040  FM=0.0 

XTEST  =XTEST+DA  LOG  (YM) 

I F ( (X  TEST*  DA  LOG  (HOG)  J.GT.  (-17  3.  ) ) F.'J=FQ  + FM*HOG 
I F ( ( XTE$T+DAL0G(H0G2  ) ) ,GT.(-173. ) )F0?=FU2+FM*H0G2 
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IF ( XTEST. OT. (-173 . ) ) FM=YN*YN1 *Z  1 

IF(  ( XTEST +DALOG(HKG)  ).GT.(-173.)  ) F = F+  F M*HKG 

IF  ( HNG.  S3 . ( 0. 0 ) ) G(J  TO  2080 

I F ( (XTEST+DALOGThNC,)  ).GT.  (-173.  ) ) F 2 = r 2 ♦FM*HNG 

-*PQN=1  .2041 2*IABS ( I XTUV(  F2,I ,7 )-IXTJV(FM,l ,7  ) - 1 XTU  V ( HNG  , 1 ,7  )+64  ) 
IFIFXPON.LT. ( 16. ) )G0  TO  2030 

KUUNT  =KOUNT  + 1 

WRITE(6,2009)cXP0\ 

2060  CONTINUE 

2090  E 1 1=0 .0 

PERCNT  = <QJNT*100. /M 

I F ( (OALOG(IB)  +DALOG  ( = ) ) T . (-173  . ) ) - 11  = 1 'A*? 

E 1 1 .?  =0  . o 

I F ( (2.*DAL0G(  I*)+PALOC ( r 2 ))  .GT.  (-173.  ) )r  1 12  = 13*1  3*F_’ 

M2 ( I * J ) =F7 

M3(I,J)=F02 

Tc(  (2.-0AL JG(^ J) ) .07. (-172. ))  -0( I , J)=M3( !, J)-F  »**“ 

MEAN=Mt AN  + E 1 1 

VAR=VAR+C I 12 

I F ( ( 7.  *OALOG(  Ell)  ) .0T.(-f  7j /)  )VAR  = /AR-m**7 
WRITE  (6,1  1 1 » ME  AN  * 7 VR,.>L<C,\T 

12  FORM A f ( • V, 6 ( 1 X » 1 POl 3 . 3 ) » 1 X , 5 % . 3 ) 

90  CONTINUE 

100  CONTINUE 

STDV=NS  T (VAR) 

WRITS (6, 1"!  ) AM, ME  AN , CT  i )V 

101  FORMAT  (•  AM  = * , 1.P01  i.6,  * * AcN  = * , 1 P01  1 . = , » IT.  .v-  7.=*,  1 ?2  1 3 .6  ) 

110  CONTINUE 

CALL  EXIT 
p'NO 

IEB318I  HIGHEST  CONDITION  UP-  «'A  n O o or . rM » 1 1 0 
IEB919I  END  OF  Jno  I'hUFOTr. 
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